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PREFACE 


This report documents the technical results obtained from navigation system 
definition studies performed under Tasks 5 and 15 of Contract NAS5-26546, 
"Tracking Data Acquisition System (TDAS) Study". Task 5 (TDAS System Archi- 
tecture) and Task 15 (Alternative User Navigation Techniques with TDAS) were 
part of a two year pre-Phase A concept definition study for TDAS as the pro- 
posed successor to TDRSS, currently under development. 

SCOPE OF WORK 


The TDAS study covers a fifteen year planning period, (1995-2005). Potential 
missions to be flown In this time frame Include free flyers, support vehicles 
and a space station/platform. Much of the TDAS requirement will be to support 
low earth orbit (LEO) missions In terms of communications, navigation and TT&C. 
Additional requirements could stem from user mission activities In higher 
(e.g., synchronous) orbits, and In support of Inter-orbital transfers of mate- 
rials and men for maintenance and repair in space, or for retrieval of plat- 
forms and experiments. 


This report Is Volume VI of nine volumes constituting th'e final report for 
the TDAS pre-Phase A study. Volume titles, largely self-explanatory, are 
given below: 


e Volume I 

• Volume II 

• Volume III 

• Volume IV 

• Vol ume V 

• Volume VI 

• Volume VII 

• Volume VIII 

• Volume IX 


Executive Summary 

TDAS User Community Characteristics 
TDAS Communications Mission Model 
TDAS Space Segment Architecture 

TDAS Ground Segment Architecture and Operations Concept 

TDAS Navigation System Architecture 

TDAS Space Technology Assessment 

TDAS Frequency Planning 

TDAS Cost Summaries 
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OVERVIEW OF VOLUME VI 


The navigation architecture aspect of this study involves examination of 
TDAS-based tracking alternatives for providing user orbit and time deter- 
mination (OD/TD). Two-way range and doppler tracking as implemented in 
TDRSS is also an alternative for ground-based navigation support with TDAS. 
However, the primary focus here will be on one-way range and doppler tracking 
methods, specifically: 

• Forward Link Beacon Tracking (FLBT) - with on-board processing of 
independent navigation signal transmissions broadcast continuously 
by TDAS spacecraft, 

• Forward Link Scheduled Tracking (FLST) - with on-board processing 
of navigation data received during scheduled TDAS forward link ser- 
vice intervals, and 

e Return Link Scheduled Tracking (RLST) - with ground-based processing 
of user-generated navigation data during scheduled TDAS return link 
service Intervals. 

This study addresses system configurations and requirements to support each 
method and assesses the potential navigation performance as a function of 
user orbit, TDAS constellation options and other parameters. Results are 
then compared with accuracy requirements in the TDAS mission model. Impacts 
of the above alternatives on both TDAS and users are evaluated and key issues/ 
tradeoffs are identified. The study also considers TDAS satellite tracking 
for two options (BRTS and VLB I ) * to identify system configuration impacts of 
the TDAS constellation options and compare potential tracking accuracy per- 
formance. 


* BRTS - Bilateration Ranging Transponder System 
VLBI - Very Long Baseline Interferometry. 
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MAJOR FINDINGS 


Results of the study are summarized In Section 6. The major conclusions 
are: 

• TDAS 8eacon Tracking (Fl.BT) will satisfy all users In the TDAS 
mission model with position accuracy requirements down to 10 m. 

• Scheduled tracking alternatives (FLST, RLST) can also meet the 
accuracy requirements except at low altitudes where performance 
Is sensitive to: 

- Drag Uncertainty 

- Frequency of Tracking and/or 

- Frequency of Navigation Data Uploads (RLST/only) 

• A two or three satellite TDAS constellation impacts performance 
as follows: 

- Selecting two satellites leads to tradeoff between coverage 
and accuracy. Increased satellite spacing improves coverage, 
but a point is reached where performance in high inclination 
orbits begins to degrade (130° spacing appears better than 
162°), 

- Selecting three satellites provides full coverage and up to 
a 2:1 advantage in navigation accuracy over two satellites. 

• Projected TDAS tracking accuracy requirements (25 m-position and 
2.5 mm/sec-velocity) can be met with VLBI tracking but not with 
a minimal* BRTS configuration. 


2 Bilateration Ranging Transponder System (BRTS) sites per TDAS spacecraft. 
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SECTION 1 
INTRODUCTION* 


1 . 1 BACKGROUND 

By the 1990's, NASA missions will require extending TDRSS capabilities for 
user tracking and orbit determination In terms of greater speed, accuracy 
and throughput efficiency. Requirements are also expected to Increase for 
on-board, real-time navigation data (orbit, time and attitude) to support 
mission operations, e.g,, data annotation, antenna pointing, rendezvous 
etc. 

A review of technology trends has identified various techniques (existing 
or under investigation) which could potentially support future orbit and/ 
or time determination needs using ground-based or quasi- autonomous methods. 
These can be broadly categorized according to the measurement technology 
employed: 

• RF Signal Detection (range, range-rate, range-difference) 

• LASER Reflection (range) 

• Optical (angle) 

A summary description of the various techniques is included In Volume II 
of the TDAS Study Report [1], 

With TDRSS, RF methods are the basis for providing user tracking services. 
The focus In this volume will be on extensions to these techniques to 
support TDAS user navigation and TDAS satellite tracking. 

The primary technique for TDRS and user orbit determination (OD) via TDRSS 
utilizes ground-derived two-way range and/or doppler data and requires a 


A glossary of acronyms is provided in Appendix F. 
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coherent forward and return link during scheduled tracking intervals (see 
[2] and Figure 1-la). This technique can also support user time determina- 
tion (TD), since clock calibration parameters (bias* drift) may be estima- 
ted simultaneously in the OD process.* An alternate tracking technique is 
also offered, which employs one-way doppler data measured at the ground 
from user* transmissions during scheduled MA return link service (see [2] 
and Figure 1-lb). User clock calibration cannot be performed with this 
method (only oscillator calibration). An advantage, however, is that re- 
turn links are more plentiful and easier to schedule than a coherent two- 
way link (20 vs 2). 

The improved coverage and faster, centralized data collection available 
with TDRSS should substantially increase the throughput efficiency, of 
ground-based OD/TD with these techniques. However, achieving the capacity 
for multi-user support with rapid turnaround or near real-time requirements 
will depend on computational enhancements at the TDRSS/OSCF. A current 
goal with hardware upgrades in process and proposed software developments 
is to achieve, by the late 1980s, a 10 minute turnaround between the end 
of a tracking pass and an updated orbit computation. [4] 

An on-board OD/TD capability would provide timely navigation data and also 
relieve TDRSS ground requirements for routine mission support. A one-way 
technique has been developed in which a TDRSS user would extract doppler data 
during scheduled MA forward link service for on-board OD only (see [5,6] 
and Figure 1-lc). Ground tests are planned for 1984 to demonstrate the 
approach and associated user equipment enhancements: doppler extractor, 

stable oscillator and navigation computer (hardware/software). [7] Further 
developments will depend on test results and user initiatives. 


Clock calibration is also possible with separate two-way measurements 
once the orbit is determined. [3] 


FIGURE 1-1: TDRSS OPTIONS FOR USER NAVIGATION 
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Wfth respect to future navigation performance, inquiries were made, as part 
of the TDAS study [lb]*, to determine potential user position and time accu- 
racy requirements. Aside from Topex-type missions*, the most stringent posi- 
tion accuracy requirement Is in the 10 m (la) range. This pertains to advanced 
resource observation type missions and a proposed space station. The most 
stringent time accuracy requirement (relative to UTC) is 1 usee (la), which 
applies to several classes of missions. 

1.2 tdas navigation architecture goals 

In accordance with the preceding discussion, the goals of the navigation 
architecture study are as follows: 

• Reduce ground requirements for routine two-way tracking support 

• Support user on-board orbit and time determination 

• Provide users timely access to ground-derived navigation data, 
and , 

9 Meet navigation accuracy requirements (> 10 m) of TDAS mission 
model . 

1.3 TASK ASSIGNMENT 

The TDRSS two-way tracking technique is a candidate for ground-based sup- 
port of user navigation in the TDAS era. This study will examine alterna- 
tive one-way tracking techniques that also support user orbit and time 
determination with TDAS. The specific alternatives to be studied are: 


Appendix A provides a summary of the TDAS mission model and navigation 
requirements. 
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• Forward Link Beacon Tracking (FLBT) - based on low power signals 
continuously broadcast by TDAS satellites to permit range and 
doppler measurements for on-board use. 

• Forward Link Scheduled Tracking (FLST) - based on ground origi- 
nated signals during scheduled contact periods which permit 
range and doppler measurements for on-board use. 

• Return Link Scheduled Tracking (RLST) - based on user generated 
signals during scheduled contact periods for ground-based range 
and doppler measurements. 

The study shall include a system level definition and requirements assess- 
ment for each configuration. An accuracy analysis based on applicable 
error models shall be performed for comparison of navigation performance 
among alternatives and against requirements in the TDAS mission model. 
Impacts on both TDAS and users and key tradeoffs will be identified and 
compared. 


1.4 OVERVIEW OF METHODOLOGY 

The methodology followed in pursuing this task is summarized in Figure 1-2. 
The first step was to define the tracking configuration in terms of functions 
and system elements needed to support each alternative. The second step 
was to consider signal processing aspects beginning with a definition of 
the tracking signal in each case, then a link performance analysis to iden- 
tify EIRP and/or G/T requirements and finally estimation of the measurement 
accuracy for range and doppler- (range-rate) data types. 

The third step was to address TDAS satellite tracking to identify possible 
system configuration impacts of TDAS constellation alternatives and to de- 
termine representative estimates of potential ephemeris accuracy. The 
analysis focused on two options: 
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FIGURE 1-2: USER ONE-WAY NAVIGATION ALTERNATIVES STIinY 
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• BRTS Tracking - the primary technique for TORS tracking [8], and 

• VLBI Tracking - proposed as a BRTS enhancement or eventual 
replacement. [9] 

The fourth step was to evaluate potential user navigation performance In 
terms of OD/TD accuracy for the three tracking alternatives. Cases for 
evaluation were derlned to assess both two and three satellite TDAS con- 
stellations, several classes of user orbits, different tracking schedules 
and two algorithms for tracking data processing (sequential and sliding 
batch) . 

The last step was to summarize the study data, compare navigation perform- 
ance results with requirements defined In the TDAS mission modal and then 
integrate the key findings In terms of a proposed system architecture to 
support user navigation functions. 

1.5 OVERVIEW OF THE REPORT 

Results of the efforts outlined above are presented in succeeding sections 
and in serveral appendices which contain much of the supporting detail. 
Section 2 describes the functions and system elements defining the config- 
uration for each tracking alternative. Section 3 and Appendix B cover the 
signal processing analysis aspects. Section 4 and Appendix C consider the 
TDAS satellite tracking analysis. Section 5 and Appendices D and E present 
results of the navigation performance evaluation. Section 6 contains the 
summary data and study conclusions. 
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SECTION 2 


USER TRACKING CONFIGURATION DEFINITION* 


The proposed TDAS architecture design [la,d] extends TDRSS data relay capa- 
bilities with enhanced SMA services for multiple access users and more SA 
channels (KSA + WSA + LSA) for single access users. This section defines 
the system configurations for supporting one-way user tracking (FLBT, FLST 
or RLST) assuming {for description purposes) that It Is provided at S-Band. 

In practice, however, the scheduled tracking alternatives (FLST or RLST) 
could also be accomodated via the SA services, as necessary. On the other 
hand, beacon tracking (FLBT) must be a generally available broadcast ser- 
vice and thus Is not compatible with the SA services. While FLBT operation 
in a totally Independent band is a possible option, frequency spectrum and 
user/system equipment considerations tend to favor S-Band reuse. 

The following two subsections present the system functions and major ele- 
ments needed for supporting each alternative. A concluding subsection 
compares tracking signals and data handling Interface requirements for 
the user, space and ground segments. 

2.1 FORWARD LINK TRACKING (FLBT/FLST) 

An overview of the major functions involved with the two forward link track- 
ing alternatives is shown In Figure 2-1. Tracking signals are uplinked from 
the ground at K^-band, relayed by one or more TDAS satellites and received 
at S-band by the user for on-board processing. The fundamental differences 
between beacon and scheduled tracking relate primarily to the signal avail- 
ability and format {discussed in Section 3). 


A glossary of acronyms Is given In Appendix F. 
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FIGURE 2-1: FORWARD LINK BEACON £ SCHEDULED TRACKING (FLBT & FLST) 

(FUNCTIONAL OVERVIEW) 





The beacon signal Is assumed to be generated by a TDAS ground terminal for 
broadcast by each TDAS satellite using a single element of the 60 element 
SMA antenna array. Users would receive continuous tracking signals while 
within a beacon antenna's ±13° field of view. The corresponding upper 
limit on user altitude for 100% coverage exceeds 3100 km for all TDAS con- 
stellations considered.* 

In the scheduled mode, signals for tracking are available only during an 
allocated contact period as part of normal SMA service. Since each TDAS 
can support two SMA forward channels [Id], up to four simultaneous users 
theoretically could be supported with a two satellite TDAS constellation 
depending on channel scheduling policy. 

A block diagram of essential elements in the forward link tracking configu- 
ration is shown in Figure 2-2. For FLBT support, a separate (beacon) sig- 
nal generator Is employed at the TDAS ground station and a dedicated K u ~ 
to-S Band repeater Is needed in each TDAS satellite [Id]. For either FLBT 
or FLST, the user requires a stable frequency standard, range and doppler 
extractors, a decoder for (ancillary) navigation data and a computing 
facllfty for on-board 0D/TD processing. Also, the required ancillary data 
(time reference information, TDAS ephemerldes, link calibrations, etc.) 

Is generated by the ground segment and incorporated in the tracking signal. 

2.2 RETURN LINK TRACKING (RLST) 

An overview of the major functions involved with this alternative is shown in 
Figure 2-3. During scheduled return link intervals, user S-band transmissions 
for tracking are relayed to the ground (at K u -band) by the assigned TDAS satel- 
lite. Range and range-rate (doppler) data are measured at the ground and ancil- 
lary data (timing**, position estimate, etc.) is decoded. The tracking data Is 
processed to determine the user's orbit and clock parameters (e.g., bias, drift). 

* Lower altitude coverage is governed by the zone of exclusion size for a 
given constellation (See [Id] for details). 

* This includes the user's current time and a sync word for range ambiguity 
resolution at the ground. 
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To support users with on-hoard navigation data requirements, the ground 
also generates orbit and time prediction data for subsequent upload to the 
user during a scheduled forward link service Interval, Such prediction 
data could comprise algorithm coefficients for Keplerian arc and/or poly- 
nominal curve fitting models [10,11] and the applicable time Intervals. 

A user receives the ground-derived navigation data and computes orbit and 
time prediction between uploads. The Interval between uploads Is an Impor- 
tant factor considered In the navigation performance evaluation (see Sec- 
tion 5), 

A block diagram of the essential elements In the return link tracking con- 
figuration Is shown In Figure 2-4. For the ground segment, the same types 
of equipment used for two-way tracking support and user orbit verification 
could also be used for RLST functions. A ground computing capability would 
also be necessary for user OD/TD and prediction data generation with ade- 
quate turnaround. The user requires a stable frequency standard, means 
for ancillary data generation and formatting, and a computing facility for 
navigation data interpretation and propagation between uploads. There 
would be no impact on the TDAS satellite configuration with this option. 

2.3 TRACKING SIGNAL AND DATA HANDLING INTERFACES 

The multiple beam antenna and switch enhancements for TDAS spacecraft pro- 
vide the capability for simultaneous, direct transmissions between the 
space segment and several ground stations [ld,e]. This provides possibil- 
ities for direct control of user spacecraft by the mission control centers 
(MCC) Instead of Interfacing through the Network Control Center (NCC) and 
White Sands (WSN) terminal as In TDRSS (see Figure 1-1). 

Figure 2-5 Illustrates options for tracking signal and navigation data 
flow with each of the one-way alternatives. Since the beacon signal for 
FLBT is a general resource, It Is assumed to originate at WSN, the assumed 
control point for TDAS spacecraft [le]. Navigation data computed on-board 
can be received by a TDAS ground terminal at the NCC directly and by MCCs 


FIGURE 2-4: RETURN LINK SCHEDULED TRACKING SYSTEM CONFIGURATION 
























with direct space/ground access. Additional Interfacility transfer of 
data can occur to support MCCs without direct space/ground access or for 
general coordination and/or verification functions. The latter Includes 
WSN which has control of TDAS spacecraft facilities (e.g., SA antenna/tele- 
scope pointing). 

With FLST (see Figure 2-5b) the user tracking signal Is Imbedded In the 
normal uplink data communication traffic so It can emanate from either the 
NCC or a cognizant MCC. Navigation data computed on-board can be distri- 
buted In the same manner discussed above for FLBT. 

With RLST (see Figure 2-5c) the user tracking signal Is Imbedded In the 
normal downlink data communication traffic, so it can be received by either 
the NCC or cognizant MCC. Ground processing for user OD/TD can occur at 
the Orbit Support Computer Facility (OSCF) or at the MCC with subsequent 
interfacility data transfer as noted above. Since the two-way tracking 
with ground-based processing is also a TDAS alternative, Figure 2 6 Is 
Included for comparison with the RLST data flow. 
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SECTION 3 


TRACKING SIGNAL CHARACTERISTICS 

The one-way alternatives for user navigation are distinguished by tracking 
signal generation, transmission, and/or processing aspects. This section 
presents: a tracking signal definition for each case (Including a candidate 

beacon signal); results of link performance analyses to assess EIRP and G/T 
tradeoffs; and estimates of the measurement accuracy for range and range- 
rate (doppler) tracking data. 

3.1 FORWARD LINK TRACKING (FLBT/FLST) 

3.1.1 Tracking Signal Definition 

On-board navigation by user spacecraft, as defined, requires that the 
received forward link signals provide the means to make range and range- 
rate measurements. Table 3-1 details the possible signal characteristics 
for two different modes: Forward Link Beacon Tracking (FLBT) and Forward 

Link Scheduled Tracking (FLST). Users may perform doppler only tracking 
of the signals to obtain range-rate values, but the signals are modulated 
by pseudonoise (PN) codes to also allow estimation of range delay. A chip 
rate of 3.0778 Mcps Is assumed here in keeping with current TDRSS practice, 
although higher chip rates are perhaps possible. The employed PN codes are 
unique to each user in the case of FLST and unique to each TDAS satellite in 
the case of FLBT. The signals are also modulated by data, providing the user 
with, as a minimum, the TDAS ephemeris and the timing information necessary 
for clock synchronization purposes. It is assumed that, in addition to on- 
board navigation capabilities, user spacecraft will also possess Viterbi 
decoders to permit convolutional coding of the forward link data and a con- 
sequent 5 dB coding gain. Current technology has already placed such decoders 
on VLSI chips. 

Forward Link Beacon Tracking : The ground-generated beacon signal is con- 

ceived to be transmitted to user spacecraft by a single element of the 
60-element TDAS SMA antenna array. Due to the large beamwidth of a single 
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element, the beacon Is available to all users within the TDAS field of view 
(±13°). In 'supporting more than one user simultaneously, however, adjust- 
ment of the transmit frequency to compensate for doppler shifting between 
TDAS and the user is not feasible. The consequent search in frequency 
required for a user to acquire the beacon signal may be significantly 
reduced by orbit predictions from prior tracking intervals. Coordination 
with other modes of tracking would provide a basis for such knowledge. 

Beacon signal acquisition is discussed further in Section 3.1.3, 

The beacon signal is assumed to comprise PN modulated BPSK data transmitted 
at a data rate of approximately 100 bps. Such a low data rate is consis- 
tent with: the EIRP provided from a single TDAS SMA antenna element, low 

gain user antennas, and system noise temperatures. The data would be assem- 
bled into frames of roughly 1000 bits, each frame containing, for example: 

• TDAS ephemeris parameters 

- Keplerian arc parameters 

- ephemeris reference time 

- age of data 

• TDAS timing and link calibration parameters 

- PN code epoch 

- polynomial correction coefficients 

- age of data 

• Frame overhead 

- frame synchronization word 

- error control field 

• Miscellaneous data 

- almanac parameters for other TDAS satellites 

- health of all TDAS satellites. 
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The data is assumed to be synchronous with the PN code, allowing bit syn- 
chronization of the data to resolve range ambiguity, e.g., [12,13]. 

Forward Link Scheduled Tracking : This mode is discussed based on data for 

navigation derived from scheduled forward link SMA transmissions. The sig- 
nal structure is assumed to be the same as in TDRSS, Separate, simultaneous 
command and range channel signals are transmitted using staggered quadriphase 
shift keying ( SQPSK) as the modulation technique. The command channel con- 
tains all data and is asynchronously modulated by a short (1023-chip) PN 
code, thus permitting rapid acquisition; the range channel is acquired sep- 
arately and is modulated by a long (256 x 1023-chip) PN code to allow the 
resolution of range ambiguity. Since data is contained only in the command 
channel, unbalanced SQPSK modulation is employed, maintaining the ratio of 
power in the command channel to that in the range channel at 10 dB. The 
command and range channel PN codes are time-synchronized so that acquisition 
of the command channel's short code assists in acquisition of the range chan- 
nel's long code. Acquisition of both the carrier frequency and PN codes Is 
further aided by doppler compensation so that the carrier frequency received 
by the user is within a predictable tolerance. 

3.1.2 Link Performance 


Assessment of the error contribution in range and range-rate measurements 
due to thermal noise depends on the received C/N 0 , the ratio of signal power 
to noise spectral density. Table 3-2 shows the TDAS SMA forward link budget 
used in this study. The TDAS satellites act as relays from the ground to 
the user spacecraft, comprising a two-link communication system. Since the 
weaker of the two links constrains the overall communication performance, 
only the satellite-to-user link is considered here. Reasonable values are 
assumed for the polarization and miscellaneous system losses; the path loss 
is calculated assuming a satellite-to-user distance of approximately 41,000 
km. The ratio C/N Q for the forward link is thus parameterized by the user 
spacecraft's antenna gain and system noise temperature (G/T) and the EIRP 
emitted by the TDAS satellite. 
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The maximum data rate that can be supported at a given bit error rate Is 
also a function of the received C/N 0 . If a desired 10"® bit error rate, 
a 5 dB gain obtained from convolutional coding, and a 3 dB margin are 
assumed, Figure 3-1 shows the achievable data rates as a function of user 
G/T, parameterized by several representative TDAS power levels. For the 
beacon case, decoding 100 bps ancillary data Implies a user G/T requirement 
> -32 dB/°K. This Is not a very stringent requirement considering that a 
0 dB antenna gain and a 1000° K receiver system noise temperature Is easily 
attainable. 

3.1.3 Beacon Signal Acquisition 

The requirements for user acquisition of the beacon signal are affected 
by Its frequency uncertainty due to doppler shifting. As an example, the 
Motorola 2nd generation TDRSS user transponder design [14] specifies a 
minimum C/N 0 of 33 dB-Hz for acquisition In 20 seconds with a probability 
of .9, assuming a frequency uncertainty of 1.5 kHz and a two channel search 
In time and frequency. Figure 3-2 shows the user G/T and TDAS EIRP required 
to surpass this threshold. With a single-element EIRP of 25.9 dBW for the 
beacon signal, the user G/T must be greater than approximately -27 dB/°K 
to meet the conditions of the Motorola design. The uncertainty In frequency 
and time for the beacon will lead to longer acquisition times, but this 
may be acceptable In the case of beacon acquisition. 

To ameliorate the signal acquisition task, future user transponder designs 
might use more than two channels to effect the PN code acquisition. Figure 
3-3 shows a 3-channel configuration. All three channels could be used In 
acquiring the beacon with a resultant decrease In acquisition time. When 
scheduled service begins, one channel may remain dedicated to tracking the 
beacon while the remaining two channels acquire the doppler compensated for- 
ward link scheduled service. Use of the three channels could be directed 
by the user's on-board computer In various handover scenarios, as suggested 
in Figure 3-4. Furthermore, the acquisition search in time and frequency 
could take advantage of past tracking data to optimally conduct the search. 
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FIGURE 3-2: TDAS SMA FORWARD LINK: PN CODE ACQUISITION THRESHOLD 
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3 , 1,4 Tracking Data Measurement Accuracy 
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From the analyses of Appendix B, measurement errors for both FLBT and FLST 
modes were estimated and are shown In Table 3-3. The thermal noise calcu- 
lations assume that the user G/T equals -27 dB/°K as a baseline value; 
code and carrier loop parameters are assumed to be In keeping with current 
design practice. It can be seen that the range error Is dominated by sys- 
tematic delay errors encountered In the ground-to-TDAS link. These are 
assumed to be roughly the same as specified for TDRSS, with a total sys- 
tematic delay of 10 meters conjectured if a routine and thorough system 
calibration capability is Implemented [15], 

Range-rate error for FLBT Is dominated by thermal noise, reflecting the 
relatively low power of the beacon signal. Increased beacon power and, of 
course, higher user G/T values would diminish the error. For FLST, phase 
noise Introduced by the ground frequency standard and the voltage-controlled 
oscillators (VCO’s) on-board the TDAS satellite are major components of the 
total range-rate error. The estimated phase noise contribution is derived 
from consideration of possible TDAS mixing frequency schemes, as detailed 
In Appendix B. The alternative of a K a -band ground-to-TDAS link versus the 
K u -band link used in TDRSS introduces additional phase noise Into range-rate 
measurements due to the larger frequency multiplies Involved, thus domina- 
ting the total range-rate error. Use of a crystal oscillator disciplined 
to the ground cesium frequency standard would, however, Improve the ground 
contribution to phase noise. To estimate the Impact, the analysis of 
Appendix B.3 was applied using the crystal oscillator phase noise spectral 
density for frequencies less than 1000 Hz and that of the cesium oscillator 
for frequencies greater than 1000 Hz as an approximation to the possible 
ground phase noise spectral density. This leads to a ground contribution 
of roughly .3 mm/s with a Ku-band pilot tone and roughly .6 mm/s with a 
Ka-band pilot tone. In comparison to the pessimistic values cited In 
Table 3-3 (1.8 and 2.0, respectively), the improvement Is significant. 
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3.2 


RETURN LINK TRACKING (RLST) 


3.2.1 Tracking Signal Definition 

As a one-way navigation option, RLST involves the ground tracking of non- 
coherent transmissions by user spacecraft using the S-band multiple access 
return service. The signal structure is assumed to be consistent with 
TDRSS DG1 Mode 2 [2] and is shown In Table 3-3. As with forward link 
tracking, an ultra-stable oscillator is required on-board the user space- 
craft to provide an accurate time and frequency standard. The user-unique 
2047-chip Gold codes on the I and Q channels are not long enough to resolve 
range ambiguity; the data must therefore be synchronously modulated by the 
PN codes to allow bit and frame synchronization at the ground receiver to 
el 1 ml Mate the ambiguity. To support one-way navigation, the data must 
contain the user's PN code epoch and timing information. 

3.2.2 Link Performance 


In order to evaluate the measurement errors in Return Link Scheduled 
Tracking due to thermal noise, the link budget of Figure 3-5 was used to 
caluclate C/N 0 . As in the calculations for forward link tracking, the 
stronger TDAS-to-ground link is neglected here. The budget assumes a G/T 
figure of -14.1 dB/°K per element of the 60-element TDAS SMA antenna array 
with the 60 elements providing a theoretical combiner gain of 17.7 dB. 

The received C/N 0 at the TDAS satellite (approximately equal to that 
received at the ground) is thus a function of the user EIRP. If a 10“5 
bit error rate, a 5 dB convolutional coding gain and a 3 dB margin are 
assumed, Figure 3-5 shows the return link achievable data rate plotted 
against user EIRP. For example, a 1000 kbps user would require an EIRP 
= 0 dBW. 


3.2.3 Tracking Data Measurement Accuracy 

Table 3-5 summarizes the RLST measurement errors estimated using the tech- 
niques of Appendix B. Salient performance characteristics of the Harris 
wide dynamics demodulator were used to model the ground receiver's response 
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to thermal noise. As a baseline value, the user spacecraft Is assumed to 
transmit an EIRP of 2 dBW, thus yielding a nominal C/N 0 of approximately 
40.4 dBW at the demodulator. Systematic errors In range and range-rate 
measurements are estimated as for the forward link tracking modes. Again, 
TDAS systematic delays dominate the range error value. 

Range-rate error Is seen to depend most on phase noise Introduced by the 
ground frequency standard. This result, however, pessimistically disregards 
the possible use in TDAS of a phase-locked crystal oscillator disciplined 
to the cesium frequency standard. If such a method were employed, the 
phase noise spectral density of the ground standard would be Improved, 
reducing Its phase noise contribution to range-rate errors. Calculations 
using the phase noise spectral density mentioned In Section 3.1.4 yield a 
ground contribution of .9 mm/s for a system using a K u -band pilot tone and 
1.5 mm/s with a K a -band pilot tone. This compares to the quoted results 
of 5.7 mm/s and 8.6 mm/s, respectively. With this approximation, thermal 
noise becomes the dominate range-rate error source for RLST. 
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SECTION 4 


TDAS TRACKING ANALYSIS 


Accurate knowledge of TDAS satellite ephemerldes Is necessary for user 
navigation support and TDAS orbit maintenance functions. This section 
presents a preliminary analysis of TDAS satellite tracking to: 

• Identify potential system configuration and operational Impacts 
of TDAS constellation alternatives, and 

• Estimate potential tracking accuracy based on BRTS and VLBI tech- 
niques. 

An overview of the elements considered in pursuing this analysis is given . 
In Figure 4-1. 

4.1 TRACKING CONFIGURATIONS 

Several constellatlon/network options Identified in the TDAS study [ld,e], 
are Illustrated in Figure 4-2. Satellite spaclngs as well as frontside/ 
backside deployments were chosen to assess possible coverage improvements 
and data distribution alternatives for servicing multiple earth stations 
in CONUS. Distributed, direct access links to/from CONUS are achieved by 
satellite crosslinks and/or multi-beam antenna enhancements on the TDAS 
satel 1 Ites. 

In terms of TDAS tracking, however, the various constellation options can 
be reduced to two situations: frontside and backside satellites. The 

following subsections discuss BRTS and VLBI tracking configurations to 
support either case. 
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4.1.1 


BRTS Tracking 


Analogous to the TDRSS implementation, two-way range and range-rate tracking 
data (R,R) are assumed to be acquired with Bilateration Ranging Transponders 
(BRT) operating at pairs of automated ground stations (one pair per TDAS). 
Tracking transmissions originate at White Sands (WSN) with replys from a 
co-located BRT and one at a geographical ly displaced site. A typical data 
rate would be 1-5 observations per minute for 5 minutes every hour, as 
expected for TDRSS. [8], 

Figure 4-3 lists some geometrically compatible station pairs* that cover 
four of the TDAS locations considered, (see also Figure 4-4) Two satel- 
lites are at TDRS locations and a third (at 100°W) is In full view of CONUS. 
The fourth is a backside satellite (at 98°E) which requires a satellite- 
satellite crosslink In the BRTS configuration if WHS is the base location 
for originating transmissions. This Impacts tracking accuracy, since meas- 
urements are also affected by uncertainties In .the frontside satellite orbit. 

The last two BRTS configurations suggested for backside TDAS tracking in 
Figure 4-3 employ two backside stations to offset the crosslink impact. In 
the first approach, one station is configured to emulate WHS for tracking 
data generation and measurement, but this would add a significant overseas 
hardware and maintenance requirements. In the second approach, each two-way 
measurement made between WHS and a backside station pair (S 2 » S 2 n ) is 
differenced with a corresponding measurement (My) made between WHS and 
the backside TDAS. The effect is to cancel systematic uncertainties In the 
crosslink and frontside space/ground link. However, this also increases 
the noise error due to the additional measurement (My), although increasing 
the data rate sufficiently may permit additional data smoothing to offset 
this. 


* 


Political issues regarding potential overseas sites were not considered. 

























4.1.2 VLB! Tracking 


The VLBI technique Is under Investigation as a possible BRTS enhancement 
[9] with potential for high accuracy TDRS tracking. The fundamental VLBI 
measurement data type Is equivalent to the difference In range (aR) between 
a signal source and two receivers displaced along a known baseline. For a 
typical measurement, wideband RF signals are received at two sites, mixed 
down to baseband and recorded over a prescribed Interval. Signal phase 
coherence and time synchronization are maintained at each site with an 
ultra-stable frequency standard. Subsequent crosscorrelation of the two 
signal streams determines the relative difference In time of arrival which 
equals {after converting to distance) the geometric range difference plus 
measurement errors {random and systematic). Investigators at GSFC using 
signals sampled from celestial radio sources by several VLBI stations 
report capabilities for AR measurement precision In the 0.3-3 cm range. 

[ 20 ] 

For the TDAS application, tracking signals assumed to originate at WHS are 
relayed by the satellite constellation to automated VLBI ground stations. 
Received signals are time tagged and buffered and, on command, returned via 
TDAS to WHS for processing. A representative measurement data rate would 
be one processed pair of VLBI observations (aR) every hour per station set. 
Nominally, one set of three stations with adequate baseline geometry would 
be needed per TDAS satellite. 

Figure 4-5 lists some geometrlcally-compatible 3 station sets that cover 
the 4 TDAS satellites considered -(see also Figure 4-4). The station sets 
are grouped In terms of long and moderate baselines, the latter of interest 
In a CONUS-based VLBI network for frontside TDAS tracking. For the backside 
satellite, no significant Impact on measurement accuracy due to the cross- 
link would be expected, since the equivalent tlme-of-arrl val differencing 
in the signal processing tends to cancel common path uncertainties. 


FIGURE 4-5: VLBI CONFIGURATIONS FOR IDAS TRACKING 







Relative errors between stations (e.g., clock synchronization and baseline 
uncertainties) are a general concern, since they can have a significant 
Impact on TDAS OD accuracy. A possible alternative to relieve stringent 
calibration requirements at each station Is to periodically Include VLB I 
differences, l.e., AVLBI measurements In the tracking process. For this, 
additional VLB I observations derived from a known celestial radio source 
are subtracted from TDAS VLBI observations made for the samo station pair. 
The effect Is to cancel station time synchronization errors and reduce 
effects of station survey errors on TDAS OD accuracy. 

Of course, the A VLBI alternative would liilrcduco additional equipment com- 
plexity at each station to also receive celestial source slgnal-i. 
handling requirements would also be Increased. Equipment cf)ii$lexf'iy afai 
operational aspects are Important Issues especially fcr VLB l 

tracking since 3 remote stations per TDAS arc Involved compa rqo to l is 4 !"' 
perhaps 2) for BRTS. 

4.2 TDAS TRACKING ACCURACY 

Information about TDAS orbits Is assumed to be provided for user OD/TD on 
a recurring basis. User navigation performance Is then a function of TDAS 
orbit uncertainties In the Interval between updates {l.e., prediction Inter- 
val errors). To assess potential orbit prediction errors, an accuracy ana- 
lysis was made of TDAS tracking with the 8RTS and VLBI techniques. The 
following subsections discuss the error modelling approach and the major 
results. Additional details are given In Appendix C. 

4.2.1 Error Model! inq 

In this study it was assumed that TDAS OD Is performed for a prescribed 
epoch via batch processing of tracking observations taken over a given data 
arc (or tracking Interval). The epoch solution propagated beyond the end 
of tracking (EOT) Is the satellite orbit state in the prediction Interval. 
Orbit prediction error covariances were evaluated using the ORAN program, 
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a linear covariance analysis tool [16, 17], ORAN was configured to compute 
satellite orbit uncertainties versus time from epoch, given the following 
Input data: 

• Initial State Vector - All TDA$ were assumed to be In 5°, circular, 
geosynchronous orbits with epoch location at the nominal longitude 
indicated in Figures 4-3 and 4-5. 

• Tracking Station Locations - Stations pertaining to a given TDAS 
and tracking technique are identified in Figures 4-3, 4 and 5 
with specific locations as tabulated In Appendix C.4. 

• Tracking Schedule - BRTS observations were assumed to occur at 
1/mln for 5 mins/hr and VLBI observations at 1/mln to 1/hr*. 

• Tracking Error Model - Values are stated in Table 4-1. All local 
and dynamic errors were treated as systematic consider error 
sources In the analysis. 

BRT.i measurement errors are consistent with values used for TORS OD analyses 
[18] assuming comparable equipment and link quality. Station errors were 
Improved by 3:1 assuming better site survey capabilities in the TDAS time- 
frame. 

VLBI measurement errors were considered at two levels: a baseline model 

consistent with values used in studies of a Deep Space Network (DSN) tracking 
application to TDRSS [19] and a reduced model which reflects achieved 
capabilities with VLBI stations observing celestial radio sources [20]. 
Station errors were taken to be the same as BRTS for the baseline model and 
an order of magnitude better for the reduced model . 


OD performance was evaluated as a function of VLBI data rate. (See Appendix C.3). 
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The tropospheric and solar pressure errors reflect modelling Inaccuracies 
and are standard values used In TDRS OD studies [18,19], although some 
Improvement may be possible. The gpavitational errors are conservative 
since further Improvements are likely by the TDAS timeframe. 

4.2.2 Error Analysis Results 

TDRS studies Indicate that the choice of tracking Interval and reference 
epoch can significantly Impact both definitive and predictive OD accuracy, 

A recent study [18] of TDRS tracking via BRTS evaluated orbit position/ 
velocity uncertainties over a 24 hour prediction Interval* as a function 
of each parameter. Results show that a judicious combination can be 
selected in which the maximum uncertainty In the prediction interval Is at 
or near a minimum. 

In the TDAS study It was decided to also consider The prediction interval 
as a selectable parameter. The basic rationale is that a sequence of 
shorter prediction Intervals with better accuracy may be traded off against 
more frequent OD processing for user navigation data updating. In both 
BRTS and VL8I analyses efforts were made to identify multiple combinations 
(epoch/tracking interval/prediction interval) which, taken together, yield 
better accuracy over a 24 hour period than any one alone. 

4. 2.2.1 BRTS Tracking . TDAS position uncertainty was computed over a 24 
hour prediction interval as a function of tracking interval and epoch for 
both frontside (100°W) and backside (98°E) cases. Figure (4-6) presents 
results plotted from the end of tracking (EOT) for several tracking inter- 
vals with a common epoch (satellite local noon) assumed in each case.** 


The 24 hour interval was considered since TDRS OD is planned to be per 
formed only once per day. 

itic 

Epoch times are stated as Greenwich Mean Time. The tracking configura 
tions are based on Figure 4-3: frontside tracking - Configurations I 

and II; backside tracking (Option 1) - Configurations III and IV (T100 
means TDAS @ 100°W is the relay to/from the backside TDAS 0 98°E). 
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FIGURE 4-6: TDAS POSITION UNCERTAINTY IN PREDICTION INTERVAL 
VS, BRTS TRACKING INTERVAL 
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In the frontside results shown In Figure (4-6a) » the peak uncertainty Is 
due primarily to solar pressure modelling error, the leading error contri- 
butor, Consequently, the peak value turns out to be an oscillatory func- 
tion of epoch with a 12 hour cycle over epoch time of day. As observed in 
Appendix C however, even the lowest peak values exceed, by at least 30%, 
the position uncertainty for the 18 hour tracking case In the 12 hour 
period after EOT (see Figure 4-6a). Also, because of the 12 hour cyclic 
property, approximately the same error profile occurs with epoch placed 
at local midnight. Thus, orbit predictions from two 18 hour tracking 
intervals with epochs at noon and midnight could be concatenated to produce 
a comparatively flat prediction error profile over 24 hours. Comparison 
of the maximum prediction uncertainties and primary contributors for contin- 
uous and concatenated cases is shown in Table 4-2, 

Fc i backside case shown in Figure 4-6b, the results are also dominated 
by solar pressure effects. However, they also reflect the frontside TDAS 
orbit uncertainty which was assumed in this analysis to be constant (at 
the 100 m level) corresponding to the concatenated 12 hour prediction 
results discussed above. No significant benefit is apparent from using a 
similar approach in this case based on 36 hour tracking, since the peak 
uncertainty over 24 hours is nearly the same as over 12 hours. A compari- 
son of the peak values and primary error contributors can be made from the 
data in Table 4-2. Corresponding peak TDAS velocity errors are in the 
8 mm/sec and 19 mm/sec range for the frontside and backside tracking cases, 
respectively. 

The preliminary observation here is that the orbit prediction uncertainty 
over 24 hours for the backside TDAS (with tracking configuration Option 1*) 
is apparently more than double that for the frontside TDAS. Certainly 
further analysis is warranted however, to evaluate possible improvements 
that alternative configurations and/or processing approaches may yield, e.g., 


* Tracking configuration options are defined in Figure 4-3. 
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• Backside tracking configuration Options 2 or 3* which utilize 
another backside station. 

• Joint data processing for frontside and backside TDAS OD that 
accounts for mutual error source correlations that an Independent 
solution (assumed here) does not. 

4. 2. 2.2 VLBI Tracking . TDAS position uncertainty was computed over a 
24 hour prediction interval as a function of the VLBI tracking Interval 
and TDAS location for two error models defined earlier in Table 4-1. One 
Is a conservative (baseline) error model and the other is a more optimistic 
(reduced) error model. 

Figure 4-7 presents results for the baseline model plotted from the end of 
tracking (EOT) with a common epoch (satellite local noon) assumed in each 
case. The results in Figure 4-7b Indicate that VLBI tracking performance 
Is relatively insensitive to TDAS location, frontside or backside. 

As shown in Appendix C, solar pressure modelling error Is a major contri- 
butor for the longer tracking intervals (> 18 hours). Since this effect 
is also cyclic with epoch time of day, essentially the same error profiles 
result with epoch selected at local midnight. Thus, if the prediction 
interval segments for 18 and 36 hour tracking are concatenated appropriately, 
a significantly lower prediction uncertainty over 24 hours could be achieved 
than with any individual case alone.** 

With the shorter tracking and prediction intervals (£ 12 hours), measurement 
and local error sources (noise, bias and station survey) are dominant ele- 
ments. As illustrated in Figure 4-7a however, the prediction uncertainty 
increases significantly as the tracking interval decreases. While a slight 
Improvement in prediction performance can be achieved by also estimating 

* Tracking configuration options are defined in Figure 4-3. 

«L*L 

Results with VLBI bias estimation are given in Appendix C.3. These show 
that prediction segments for 18 and 12 hour tracking segments would give 
somewhat better accuracy that the 18 and 36 hour combination. 
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the VLBI measurement bias, no really significant benefit is apparent from 
using shorter tracking and prediction intervals, when VLBI measurement and 
local errors are at the levels assumed in the baseline model. However, 
if a VLBI implementation can be realized with these errors at the level assumed 
in the reduced model, a different situation emerges. 

Figure 4-8 compares TDAS position uncertainties for short tracking and 
prediction Intervals for both error models. In the latter case, the lowest 
position uncertainty over a given prediction Interval occurs with 6 hour 
tracking. As observed in Appendix C, this is essentially the cross-over 
region between tropospheric and solar pressure effects with the latter 
becoming significant for longer tracking/prediction Intervals. Thus, a 
sequence of short prediction intervals (e.g., 1 hour) could provide a 
position uncertainty in the range of 10 - 15m. 

Figure 4-9 shows examples of TDAS position uncertainty from using multiple 
(concatenated) prediction segments over a 24 hour period. In the baseline 
case (Figure 4-9a )* the uncertainty ranges from 50 - 80m using 4 segments 
based on 18 and 36 hour tracking or 50 - 70m based on 18 and 12 hour track- 
ing with VLBI measurement bias estimated. In the reduced case (Figure 4-9b)* 
the error ranges from 10 - 12m using 24 one hour segments based on 6 hour 
tracking Intervals with epochs spaced one hour apart, A second curve based 
on 5 % tropospheric error shows the improvement (to 5 - 8m over 24 hours) 
if better tropospheric error correction can be achieved. A comparison of 
the peak errors and primary error contributors in each case is given in 
Table 4-3. Maximum TDAS velocity errors corresponding to these position 
error results are at the 6 mm/sec and 1 mm/sec levels for the baseline 
and reduced error models, respectively. 


The assumed VLBI data rate of 1/hr (one AR pair/hr) is more representa- 
tive that the 1/min rate assumed in Figures 4-7 and 4-8. Nevertheless, 
the computed results for these tracking intervals are Identical with 
either rate. Tables C-l and C-2 in Appendix C show the impact on pre- 
diction uncertainty for other tracking Intervals. 
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TABLE 4-3 

MAXIMUM BAS POSITION ERROR IN 24 HOUR PREDICTION INTERVAL WITH VLBI TRACKING 











The foregoing results Illustrate the significant Improvement In TDAS OD 
accuracy that the VLBI technique could potentially provide, Moreover, this 
would apply to backside TDAS tracking as well, since crosslink range uncer- 
tainties are Inherently cancelled out. It must be emphasized however, that 
achieving a 10m capability Implies having measurement and station survey 
uncertainties at 4 level consistent with the reduced error model and using 
appropriate batch tracking/prediction Interval combinations {or equivalent 
scheme) for OD processing. Thus, Important areas for further study include.* 

• VLBI measurement precision (noise) with TDAS-based signals, 

• Alternatives for VLBI tracking network calibration (time synchro- 
nization and baselines) 

- via a VLB I measurements (see discussion In Section 4.1.2) or 

- via terrestrial microwave links* and a geodetic survey 
receiver (e.g. , GPS [21]) , 

• OD processing options f 

1 

„ • 

- Algorithm (batch vs sequential) , * * 

- Update rate vs prediction interval. 

m 


* This is applicable to CONUS-based tracking configurations (e.g., a con- 
nected element Interferometer network) currently under study for TDRS 
[9]. Some preliminary results for TDAS are given in Appendix C.2.3. 
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USER NAVIGATION FgRFWMANCE EVALUATION 

User navigation via TDAS Is based on computing orbits and time from range 
and range-rate measurements and known TDAS orbits. This section presents 
an evaluation of potential navigation performance In terms of OD/TD accur- 
acy for the three one-way tracking alternatives defined In Section 2. 

5.1 OVERVIEW OF APPROACH 

Figure 5-1 gives an overview of the elements considered In the analysis. 
Cases for evaluation were defined to assess several types of user orbits 
and the Impact of TDAS constellation options, different tracking schedules 
and two algorithms for tracking data processing (sequential and sliding 
batch). 

To evaluate user OD/TD accuracy two error analysis programs were employed: 

c SEA Program [221 - for tracking based on sequential processing 
of measurement data, and 

• RDGTDS [23] - for tracking based on batch processing of measure- 
ment data. 

Given a tracking schedule and nominal TDAS and user orbits each program 
computes OD/TD error covariances and sensitivities versus time with respect 
to measurement noise and other applicable error sources. The following 
subsections discuss the tracking configurations and methodology, the 
error modelling approach and the major results. Further detailed results 
are contained in Appendices D and E. 
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5.2 


TRACKING CONFIGURATIONS AND METHODOLOGY 


TDAS constellation and user orbit types considered in the performance 
evaluation are defined next. Then* tracking data ^eduling and processing 
options assumed for the three tracking alternative- IT, FLST and RLST) 
are discussed. 

5.2.1 TDAS Constellations 


The various constel latl on/network options considered in the TDAS study 
were discussed in Section 4. For purposes of this evaluation the three 
constellations shown in Figure 5-2 were selected. Option 1 Is analogous 
to TDRSS with two satellites spaced 130° apart which provides 85-100% 
coverage at altitudes down to 200 km.* Option 2 also uses two satellites, 
but with the maximum allowable spacing, 162°**, which yields 98-100% 
coverage. Option 3 has three satellites, two deployed as in Option 1 
and a third on the backside, which together provide 100% coverage. 

5.2.2 User Orbits 


Six orbits were selected to compare navigation performance for various 
‘orbit altitudes as shown in Figure 5-3. As is evident from Appendix A, 
(Table A-l), they are also indicative of the majority of potential 
user orbits in the TDAS mission model. 

The low altitude (high drag) orbit types are of interest to determine 
whether more frequent tracking data, available with FLBT, is of significant 
benefit. The high and low inclination orbit types are of interest, since 
their coverage and geometrical properties can differ significantly. 

In the high inclination case, the % coverage is generally greater, but 
opportunities for (good) doppler tracking can be reduced. This occurs 
whenever the user orbit normal points sufficiently toward a TDAS, thereby 

* A. detailed discussion of coverage characteristics as a function of 
satellite spacing is given in [Id]. 

^.||i 

This is to avoid earth occulation of the TDAS-TDAS crosslink. 
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FIGURE 5-2: TDAS CONSTELLATION OPTIONS CONSIDERED 
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reducing the range change over a pass. In the limit, there Is no range 
change (zero doppler), yet full visibility. 

5.2.3 Tracking Data Processing and Scheduling Methodology 

One-way range and range-rate (R,R) are the fundamental data types employed 
for user OD/TD with each tracking alternative. System configurations 
and technical considerations for deriving these data were covered In 
Sections 2 and 3. This discussion pertains to measurement scheduling 
and processing aspects that were assumed for the performance evaluation. 

Two algorithms for tracking data processing were considered: sequential 
and sliding batch. In the sequential case, estimated parameters (e.g., 
position & velocity; clock bias & drift) are updated after each measurement 
and propagated forward between measurements and/or tracking passes. 

For this analysis an extended Kalman filter was assumed to be employed. 

In the sliding batch case, selected parameters are estimated for a pre- 
scribed epoch via standard weighted least squares processing of tracking 
observations taken over a given data span (T). Current estimates are 
derived by propagating the epoch solution over a specified interval (P) 
beyond the end of tracking and whatever computation/data handling interval 
(C) is required. In other words, current estimates are predictions in 
the interval, C to C+P , beyond the end of tracking. For each successive 
batch the process is repeated with the epoch and tracking interval advanced 
by P.* 

Table 5-1 lists the tracking algorithm and measurement scheduling options 
considered in the performance evaluation for each alternative. Only 
sequential processing was assumed for FLBT, since the sliding batch approach 
was felt to be inappropriate for on-board OD/TD in view of the tracking 
data volume and associated computation/data handling requirements. For 
FLST and RLST, two tracking schedule types were considered based on 
acquiring data every other orbit (Schedule I) and every orbit (Schedule II). 

* Since P>C, predictions for all batches extend at least 2C beyond the 
data span. In a sense, 2C is a measure of estimate "staleness". 
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TABLE 5 1: TRACKING ALGORITHM & MEASUREMENT SCHEDULING OPTIONS CONSIDERED 







For sliding batch uses, the .tracking data span was considered as 
a parameter to be selected with possible values ranging from 6-24 hours. 

5.3 TRACKING ERROR MODELING AND COMPUTATION 

In this study user navigation was evaluated from OD/TD error covariance 
computations using the SEA and RDGTDS programs [22, 23]. Figure 5-4 
illustrates some hypothetical error profiles which Indicate the peak 
errors and time intervals of interest for the options considered: 

• Sequential Processing - For FLBT and FIST, user OD/TD accuracy 
was defined as the peak error occurlng over 24 hours with track- 
ing Intervals scheduled according to Table 5-1. For RLST, 1 

this criterion also applies to ground-based operations, but 
not necessarily for on-board navigation. Since navigation 
data updates would be available only on a recurring basis, 
user OD/TD accuracy was defined instead as the peak error 
In a prediction-only mode between uploads. 

* ! 

i 

9 Sliding Batch Processing - For FLST and RLST, user OD/TD accuracy \ ) 

was defined as the peak error occurring over 24 hours in a series 1 ' ; 

of batch-derived prediction Intervals (P), see Figure 5-4. This . i 

is defined by the number of batches/day for FLST and by the ' ' < 

I 

navigation data upload rate for RLST. In either case, the mini- 
mum interval for P Is governed by the OD/TD computation/data j 1 

handling time (C) from the end of any tracking span.* 

To obtain user OD/TD uncertainties as a function of various parameters 
the SEA program was utilized for sequential processing cases and the . 

RDGTDS program for batch processing. Error covariances and budgets were [ . < 

computed given the following input data (see Figure 5-5): 

w , * l \ 

In Figure 5-4, C=10 minutes (for RLST) is based on anticipated ground- *" ■ 

system processing capabilities [4]; C=90 minutes (for FLST) is based 
on LSI-11 type on-board computing capabilities planned for evaluation > 

in a TDRSS forward link doppler tracking experiment [7]. 
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User Location - A user was assumed to be In one of the orbits 
(A-F) Indicated In Figure 5-2 with an Initial location at 0°N, 
0®E for the orbit selected. 

TDAS Location - All TDAS were assumed to be In 5° circular, 
geosynchronous orbits with Initial location at the nominal 
longitude Indicated In Figure 5-1 for the selected constellation 
option (1, 2, or 3). 

Tracking Schedule - (R,R) measurement times were set up to 
occur at the rate defined In Table 5-1. 

Tracking Model Parameters - Values are stated In Table 5-2 


TDAS users were assumed to estimate, as a minimum, 8 basic parameters: 

3 position and 3 velocity states and two clock, states, bias and drift. 
Unestimated parameters were treated as systematic (consider) error sources 
In the analysis, Since the analysis is linear, the results for any parti- 
cular error source may be scaled up or down to note the impact of different 
a priori uncertainties. 


Errors due to atmospheric drag were modeled with two Interpretations of 
uncertainty In the drag coefficient (Cq) . For the highest altitude orbits 
Co was assumed not to be estimated, so 25% of the nominal value was taken 
as the consider error, a typical choice [24,33,34]. At altitudes where 
effects may be significant (e.g., orbits A,8,D,E) Cq was assumed to be 
estimated, but Imperfectly. The residual component was treated as a 
consider error with a la value equal to 2.5% of the nominal Cq. In other 
words reducing the a priori uncertainty by estimating Cq was assumed only 
90% effective. 

Values for uncertainties In the other dynamic parameters (GM, gravitational 
harmonics and solar pressure coefficient) are undoubtedly conservative for 
the TDAS time frame. Sequential processing results are based on the GEM9 
error model [35], resident, in the SEA program. Sliding batch results were 
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I TABLE 5-2: TRACKING MODEL PARAMETERS 



PARAMETER 

A PRIORI UNCERTAINTY 

ESTIMATED 

USER (H,C,L 
ORBIT | H,C,L 

USER (B 
CLOCK ( B 

500 M 

l m/sec 

1 MSEC 

200 NSEC/SEC 

UNESTIMATED 
/SYSTEMATIC^ 
\ ERRORS ) 

USER ORBITS 
DRAG C ' F 

25* 

lc D7 USER ORBITS 
A, B, D, E 

2.54* 

GRAV. CONST, GM 

0.25 PPM 

GRAY. HARMONICS (LUMPED) 
( 12 x 12) 

100* GEM9 ERROR OR 
• (GEM9-GEM7) 

SOLAR RADIATION, C R 

10* 

SYSTEM BIASES ? 

R 

10 M 

i mm/sec 

TDAS Jh,C,L 
ORBIT (H,C,L 

25,23,40 M 
3, 2, 3 MM/SEC 

USER OSC. B 

DRIFT B 

lo-io parts/day 

RANDOM 

MEASUREMENT 

ERRORS 

RANGE ° R 

5 M 

RANGE RATE a R 

5 MM/SEC 

PARAMETER 

VALUE 

TRACKING 

PROCESSOR 

TUNING 

PARAMETERS 

SEQUENTIAL! 

“ USER VEL. STATE NOISE 
- CLOCK RATE STATE NOISE 

ORBIT! 

(io-o m 2 /sec 3 (a,d) 
(10-12 H 2/ sec; 3 (b, C, E, F) 

10’ 6 nsec 2/ sec 3 (all) 

Mir*! 

- TRACKING INTERVAL 

6, 12, 18, 24 HOURS 


* DRAG COEFFICIENT (c D ) IS ASSUMED TO HE PARTIALLY ESTIMATED WITH A RESIDUAL 
UNCERTAINTY TREATED AS A CONSIDER PARAMETER. 

(NOMINAL Cq ■ 2.0, USER AREA/WEIGHT * .00272 mZ/KG) 
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derived with the GEM9-GEM7 error model, the closest option resident In 
ROGTDS. 

Constant values assumed for TDAS orbit uncertainty correspond to lo errors 
of approximately 50 m In position and 5 mm/sec In velocity, Results given 
In Section 4.2 Indicate that this accuracy level would still be fairly 
conservative based on VLB I tracking for TDAS. It would be optimistic by 
at least 2:1 for BRTS tracking*, 

Errors In the tracking measurements (R,R) were defined In terms of random 
errors and system biases. Values used for the random errors and the range 
bias error are representative of the results given In Section 3 (see Table 
3-3), The 1 mm/sec for range-rate bias was Included only to observe poten- 
tial sensitivity. 

Frequency drift In the user's reference oscillator appears as a doppler 

rate error which affects range-rate measurement accuracy and as a clock 

#* 

bias acceleration error (B) which affect}; ranging accuracy. Oscillator 
drift was defined as a consider error with a la value of 10**^ parts/ 
day, a level consistent with currently existing quartz oscillator tech- 
nology. [32] 

* 

Other modelling parameters "involved In the analysis pertain to tracking 
data processor tuning for achieving "best" performance In some sense. In 
the sequential case, state noise parameters are used as an artifice for 
adjusting the Kalman filter gains to control the weight given to prior 
estimates. The objective is to achieve a balance between uncertainties 
introduced by new measurements and those caused by propagating prior esti- 
mates with an imperfect dynamical model. The values shown In Table 5-2 
were, found to give good results In various test runs** and were adopted 

* This assumes the minimal BRTS tracking configuration using two stations 
per TDAS (see Figure 4-3). 

** As Implemented in the SEA program [22], state noise has no connection 
with physical disturbance phenomena acting on the user spacecraft or 
clock, although In theory, there could be. 
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as a baseline set for most cases considered, In the sliding batch case, 
the tracking Interval choice can be made with the same objective in mind, 
Four durations given In Table 5-2 were considered in each case and the one 
yielding the lowest peak error In a 6 hour prediction Interval from end of 
tracking was selected, 

5.4 NAVIGATION PERFORMANCE RESULTS 


TDAS user navigation performance was evaluated In terms of potential OD/TD 
accuracy for each of the one-way tracking alternatives. Results were 
obtained as a function of user orbit type, TDAS constellation option, 
tracking schedule and two processing algorithms. This section presents a 
summary and Interpretation of the significant data. More detailed Informa- 
tion can be found in Appendices D and E. 

5.4.1 Evaluation Cases 

Potential cases for evaluation were defined by the six user orbit types 
In Figure 5-2 and three constellation options In Figure 5-3. For each 
case, five subcases were defined based on processor type (sequential or 
sliding batch) and tracking schedule (I, II or beacon).* Altogether, 90 
combinations were identified and 48 of these were evaluated as indicated 
In Table 5-3. 


For each case, user position and time error profiles were computed for 
a 36 hour period (the last 12 hours In predict-only mode). From these 
the peak errors for assessing performance were determined as defined 
in Section 5,3 and Figure 5.4. Appendix D presents error profiles and 
corresponding tracking schedules for sequential processing. Appendix E 
has corresponding results for sliding batch processing. Results for 
all cases were based on the nominal model parameters given in Table 5-2 
except for certain processor tuning parameters noted specifically on 
error profile plots in Appendices D and E. 


Recall Table 5-1. 
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TABLE 5-3: CASES CONSIDERED FOR NAVIGATION PERFORMANCE EVALUATION 
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5.4.2 Sequential Data Processing Results 

Navigation performance data based on the error analysis results in Appendix D 
are presented In terms of user position and time accuracies as summarized 
In Figures 5-6 and 5-8. Corresponding position and time error budgets 
are given in Figures 5-7 and 5-9. 

5. 4. 2.1 User Position Accuracy . Figure 5-6 illustrates that in low 
altitude orbits where drag is a factor, performance depends heavily 
on frequent data availability. Beacon tracking Is clearly superior 
‘in this respect. At higher altitudes this is not the case and all three 
tracking alternatives can given comparable performance. 

With respect to the various TDAS constellations, Option 3 gives typically 
better performance due to better geometric distribution of the tracking 
data available with three satellites. In the two satellite constellations, 
Option 1 with (130° spacing) is better in the high inclination orbits, 
while Option 2 (with 162° spacing) is better in low inclination orbits. 
Although Option 2 provides nearly full coverage (> 98%), the performance 
for high inclination users is sensitive to poor geometry conditions 
for doppler tracking. This occurs twice a day when the user orbit normal 
points in the general direction of each TDAS. Consequently a tradeoff 
exists between maximum TDAS spacing (162*) for best coverage and lower 
spacing to achieve better 00 performance. 

figure 5-7 indicates the major error components affecting 0D performance. 
Those due to unestimated (consider) parameters are: gravitational harmonics 

(H), drag (D) and TDAS ephemerldes (E). The noise component (N) is 
the random error in user position which is controllable by filter tuning 
parameters.* Ideally the filter should be tuned to reduce the impact 

* By adjusting filter gains to partially deweigh prior estimates random 
errors arising from measurement and a priori uncertainties do not 
continuously decrease as the number of measurements increases. Rather 
a "noise floor'' (not necessarily constant over time) is reached, 
which is controllable using velocity state noise tuning parameters. 

[22, 37] 
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FIGURE 5-6: IDAS USER POSITION ACCURACY VS TRACKING ALTERNATIVES 

(SEQUENTIAL PROCESSING) 
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of the dominant consider parameter(s) without excessively Increasing 
the noise component. Figure 5-7 suggests that the modelled filter Is 
somewhat overtuned for the 28°, 200 km orbit, about right for the 97°, 

200 km orbit and undertuned for the higher altitude orbits. 

In the high Inclination, high altitude orbits (98°, 600/1000 km) gravitational 
harmonic modelling uncertainty Is typically the dominant error contributor. 
This Is also the case at low altitude with beacon tracking since filter 
tuning acts more to suppress the drag contribution. At low altitudes 
with scheduled tracking, drag is the dominant contributor In all cases. 

TDAS ephemerls error is typically a secondary contributor. All of the 
other consider error sources Identified In Table 5-2 are comparatively 
Insignificant Including user oscillator drift. 

Also Indicated in Figures 5-6 and 5-7 are user requirements extracted 
from the TDAS mission model (see Appendix A). Most would be met with 
beacon tracking and in many cases with scheduled tracking based on error 
analysis results using the assumed models. Meeting the more stringent 
requirements would Imply some improvement in key error sources. For 
example, In the 93°, 600 km orbit with TDAS Constellation Option 1, 
the impact of gravitational harmonic modelling errors would have to 
decrease by at least 5:1 and TDAS ephemerls error by 2:1. Filter re- 
tuning would decrease the noise component so the net effect of all con- 
tributors would be < 10 m. Such reductions could also lead to meeting 
the most stringent requirement (30 m) at a 'lower altitudes assuming 
no change in the modelled drag uncertainty. 

5. 4. 2. 2 User Time Accuracy 

Figure 5-8 indicates that TD performance is uniformly better with beacon 
tracking (< .25 usee) than with scheduled tracking in all four orbits. 
Nevertheless, performance with the scheduled alternatives is sufficient 
to meet the most stringent user time requirement (1 psec) Identified 
in the TDAS mission model. 
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FIGURE 5“8: TDAS USER TIME ACCURACY VS TRACKING ALTERNATIVES 

(SEQUENTIAL PROCESSING) 
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Figure 5-9 shows the error contributors affecting TO performance. With 
FBIT, TDAS ephemeris errors tend to dominate. With FLST and RLST, in 
the high altitude orbits, gravitational harmonic errors are also a primary 
contributor. In the lower altitude orbits the TD errors due to the 
drag uncertainty increase but are still secondary. User oscillator 
drift is also a secondary source In all cases except for RLST If the 
Interval between navigation data uploads begins to exceed 3 hours. 

5.4.3 Sliding 8atch Data Processing Results ’ 

Navigation performance data based on the error analysis results in Appendix E 
are presented In terms of user position and time accuracies as summarized 
in Figures 5-10 and 5-12. Corresponding position and time error budgets 
are given in Figures 5-11 and 5-12. FLST results are shown for batch 
Intervals (P) of 1.5 and 3 hours.* RLST results correspond to navigation 
data uploads at ^ 1, 2 or 3 orbit Intervals. 

5. 4. 3.1 User Position Accuracy. Figure 5-10 indicates that user position 
errors in low altitude orbits are substantially higher than the sequential 
provessing results. At higher altitudes however, 00 performance is 
comparable. In either case no significant OD performance difference 
was observed with respect to the two satellite constellation options 
(1 and 2) studied. 

Figure 5-11 indicates the major error contributors affecting 00 performance. 
At low altitudes drag is dominant even though a tracking interval of 
only 6 hours was used. For the higher altitude orbits a longer tracking 
interval (18 hours) was found to give better performance; consequently 
user oscillator drift became a significant contributor along with gravi- 
tational harmonic error. 


* See Figure 5-4. The assumed computation/data handling Interval (C) 
was 1.5 hours for FLST and 10 minutes for RLST as discussed in Section 5.3. 
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With respect to meeting user requirements these alternatives would not 
be good candidates In low altitude orbits. However, at the higher altitudes 
all but the most stringent requirement (10 m) could be met based on 
the error analysis results. Mettlng the 10M requirement would Imply 
some Improvement In key error sources, e.g., a decrease In gravitational 
harmonic errors by at least 4:1, TDAS ephemerls error by 1. 5-2:1 and 
user oscillator drift by 5:1. 

5. 4. 3. 2 User Time Accuracy. Figure 5-12 Indicates that TO performance 
Is at best 2.5-3 nsec. This Is significantly lower than with sequential 
processing and would not meet the most stringent time requirement of 
1 psec. 

Figure 5-13 shows that the dominant error contributors are the same 
as for user position error — drag at low altitudes and user oscillator 
drift at high altitudes. 

5.4,4 Observations 


A summary of user navigation performance based on the error analysis 
results is given In Table 5-4. In addition, the following observations 
can be made: 

• Sequential Processing 

- At low altitudes (e.g. £400 km) beacon tracking { FLBT ) 
provides significantly better orbit determination (OD) 
accuracy than scheduled alternatives (FLST, RLST), The 
latter are both highly sensitive to tracking frequency 
and RLST to NAV upload rate. 

- At higher altitudes all three alternatives provide compair- 
able OD performance. 
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- Performance projections based on error source Improvements 
Indicate that FLBT could satisfy all users In the TDAS 
mission model with position accuracy reqts. down to 10m. 

For this accuracy error sources requiring Improvement 
Include; 

— gravitational harmonic modeling (> 5:1 reduction 
In GEM-9 errors) 

— TDAS tracking accuracy (to 25m - pos., 2.5 mm/s - 
vel .) 

- Time determination (TD) Is better than 0.25 psec with 
FLBT and better than 1 psec with FLST and RLST. 

• Sliding Batch Processing 

- Results for FLST and RLST Indicate that compared to sequential 
processing; 

— OD performance tends to be worse for low altitude 
orbits {< 400 km) where drag Is the dominant error 
source and comparable to or slightly better In higher 
orbits where user oscillator drift Is also 4 major 
error source. 

— TD performance Is worse (> 2.5 psec) for all orbits 
considered and particularly sensitive to the tracking 
Interval duration. 
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SECTION 6 


STUDY SUMMARY AND CONCLUSIONS 


The TDAS navigation architecture study has focused on three TDAS-based 
one-way tracking alternatives for providing user orbit and time determi- 
nation (OD/TD). This section presents a functional overview and require- 
ments summary for each alternative, summaries of potential user navigation 
performance and TDAS tracking accuracy, and the conclusions reached from 
the study. 


6.1 FUNCTIONAL OVERVIEWS AND REQUIREMENTS SUMMARY { 2 , 3)* 

The proposed TDAS architecture design [la] would extend TDRSS navigation 
and data relay capabilities to support: 

• beacon tracking (FLBT) via an Independent beacon signal, and 

• scheduled one-way tracking {FLST, RLST) via enhanced multiple 
and single access (MA, SA) services. 

A functional overview of each alternative is presented in Figure 6-1. 

The beacon signal Is assumed to be generated by a TDAS ground terminal for 
broadcast by each TDAS satellite using a single element of the S-band MA 
antenna array. Users would receive continuous tracking signals while within 
a beacon antenna's ±13° field of view. The corresponding upper limit on 
user altitude for 100% coverage exceeds 3100 km for all TDAS constellations 
considered. Lower altitude coverage Is governed by the zone of exclusion 
(ZOE) size for a given constellation. 

In the scheduled modes, signals for tracking are available only during an 
allocated contact period as part of normal MA or SA service. Since each 

* Numbers in brackets ( ) refer to preceding sections with full details. 
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FIGURE 6-1: ONE-WAY NAVIGATION ALTERNATIVES - FUNCTIONAL OVERVIEW 

FORWARD LINK DEACON TRACKING (FUT) FORWARD LINK SCHEDULED TRACKING (FLST) 
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TOAS can support two S-Band MA forward channels and eight SA channels, 
up to 20 simultaneous users theoretically could be supported for FLST with 
a two satellite TDAS constellation depending on channel scheduling policy. 
Similarly, for RLST up to 20 simultaneous MA users theoretically could be 
accomodated, since each TDAS can support ten return S-Band MA channels. 

Table 6-1 lists various system requirements for supporting each 1-way 
alternative. Table 6-2 summarizes characteristics of the tracking signals 
assumed In each case Including estimates of metric tracking data (R, R) 
accuracies based on transmission link performance analyses. 

The multiple beam antenna and switch enhancements for TDAS spacecraft pro- 
vide the capability for simultaneous, direct transmissions between the 
space segment and several ground stations [ld,e]. This provides possibil- 
ities for direct control of user spacecraft by the mission control centers 
(MCC) instead of Interfacing through the White Sands (WSN) terminal and/or 
Network Control Center (NCC) as In TDRSS. 

Figure 6-2 Illustrates options for tracking signal and navigation data 
flow with each of the one-way alternatives. Since the beacon signal for 
FLBT is a general resource, it is assumed to originate at WSN, the assumed 
control point for TDAS spacecraft. Navigation data computed on-board can 
be received by a TDAS ground terminal at the NCC directly and by MCCs with 
direct space/ground access. Additional interfacility transfer of data can 
occur to support MCCs without direct space/ground access or for general 
coordination and/or verification functions. The latter Includes WSN which 
has control of TDAS spacecraft facilities (e.g., SA antenna/telescope 
pointing) . 

With FLST (see Figure 6-2b) the user tracking signal is imbedded in the 
normal upl ink data communication traffic so it can emanate from either the 
NCC or a cognizant MCC. Navigation data computed on-board can be distri- 
buted in the same manner discussed above for FLBT. 
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FIGURE 6-Zj TRACKING SIGNAL AND DATA HANDLING INTERFACES (ONE-HAY TRACKING ALTERNATIVES) 
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With RLST (see Figure 6-2c) the user tracking signal is imbedded in the 
normal downlink data communication traffic, so it can be received by either 
the NCC or cognizant MCC. Ground processing for user OD/TD can occur at 
the Orbit Support Computing Facility (OSCF) or at the MCC with subsequent 
interfacility data transfer as noted above. Since two-way tracking with 
ground-based processing is also a TOAS alternative, Figure 6-3 provides a 
comparison with the RLST data flow. 
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6.2 


USER NAVIGATION PERFORMANCE SUMMARY 


To assess the potential navigation performance with each of the one-way 
alternatives, user OD/TD accuracy was evaluated as a function of various 
parameters and compared with requirements In the TDAS mission model. 

Figure 6-4 gives an overview of the elements Involved In the analysis. 

User Orbits and TDAS Constellations (5.2) - Options considered are shown 
in Figure 6-5. Constellation Option 1 is analogous to TDRSS with two 
satellites spaced 130° apart which provides 85-100% coverage at altitudes 
down to 200 km. Option 2 also uses two satellites, but with the maximum 
allowable spacing, 162°, which yields 98-100% coverage. Option 3 has three 
satellites, two deployed as In Option 1 and a third on the backside, which 
together provide 100% coverage. 

The low altitude (high drag) user orbit types are of Interest to determine 
whether more frequent tracking data, available with FLBT, is of significant 
benefit. The high and low Inclination orbit types are of Interest, since 
their coverage and geometrical properties can differ significantly. 

Tracking Schedules/Processlng Algorithms (5.2) - For the FLBT mode, the 
schedule was assumed to be continuous for a selected TDAS with metric 
tracking data (R,R) processed sequentially. For the FLST and RLST modes, 
the tracking schedule Impact was assessed based on tracking every orbit 
and every other orbit during a 10 min pass/TDAS. Also, navigation perfor- 
mance was evaluated for both sequential and sliding batch* data processing. 

Error Modelling - User OD/TD error was computed via covariance analysis 
programs [22,23] given nominal TDAS and user orbits, a tracking schedule 
and processing approach, and appropriate models of the tracking error 
sources. The latter comprised measurement errors and various systematic 
errors, e.g., gravitational harmonic modelling error, drag modelling error, 
TDAS ephemeris error, user oscillator drift affects, etc. 

* Details are given in Section 5. 
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FIGURE 6-<li NAVIGATION PERFORMANCE EVALUATION - OVERVIEW 
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FIGURE S— 5 1 USER ORBITS AND TDAS CONSTELLATIONS CONSIDERED 
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Table 6-3 lists the baseline modelling assumptions made for these partic- 
ular error sources and their relative Impact on achieving a 40 M OD accu- 
racy with beacon tracking In the user orbits considered, in this case 
gravitational harmonic error dominates, so TDAS ephemeris error is only 
a secondary effect. The projected requirements indicate the gravitational 
model Improvement required (In terms of GEM-9 errors) to achieve 10 M and 
30 M accuracies in the noted orbits. At this level, the modeled TDAS 
ephemeris error would also be significant so a reduction of *2:1 is indi- 
cated. 

Navigation Performance With Sequential Processing (5.4.2) - Figure 6-6 
shows the OD performance (max error over 24 hours) with beacon tracking 
for each TDAS constellation option along with user requirements identified 
in the TDAS mission model (see Appendix A). Higher errors tend to occur 
for Option 2 in near polar orbits particularly at low altitudes due to 
recurring poor geometry.* 

Error analysis results correspond to the basic modelling assumptions and 
the projected results to the Improved errors Indicated in Table 6-3. Thus, 
beacon tracking is projected to satisfy all user requirements in the TDAS 
mission model down to 10 M. 

Figure 6-7 shows the navigation performance with FLST for two assumed 
tracking schedules and sequential processing. In contrast to beacon 
tracking, the performance in low orbits is drag dominated and substantially 
worse. At higher altitudes, results are comparable. 

Figure 6-8 shows the navigation performance with RLST for the same tracking 
schedules and two different nav data upload rates - 1.5 hours and 3 hours. 
Again, performance is worse in low orbits (drag dominated), but is compara- 
ble at higher altitudes. 

* This occurs about twice a day when the orbit normal points toward a TDAS 
thereby limiting the range change over the pass (poor doppler character- 
istic). 
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FIGURE 6-6 

USER NAVIGATION PERFORMANCE .SUMMARY - BEACON TRACKING 
(SEQUENTIAL PROCESSING) 
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TD accuracy results* Indices that performance Is uniformly better with 
beacon tracking (< .25 u?ec) than with scheduled tracking In all orbits 
considered. Nevertheless, performance with the scheduled alternatives Is 
sufficient to meet the most stringent user time requirement (1 usee) Iden- 
tified In the TDAS mission model. 

Navigation Performance With Sliding Batch Processing (5.4.3) - Results for 
FLST and RLST Indicate that compared to sequential processing: 

• OD performance tends to be worse for low altitude orbits 

{< 400 km) where drag Is the dominant error source and compara- 
ble to or slightly better In higher orbits where user oscillator 
drift is also a major error source. 

• TD performance Is worse (>, 2.5 usee) for all orbits considered 
and particularly sensitive to the tracking interval duration. 


* See Section 5.4. 
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6.3 


TDAS TRACKING ANALYSIS SUMMARY 


Information about TDAS orbits is assumed to be provided for user OD/TD on 
a recurring basis. User navigation performance is then a function of TDAS 
orbit uncertainties in the Interval between updates (i.e., prediction 
interval errors). To assess potential orbit prediction errors, an accu- 
racy analysis was made of TDAS tracking with the BRTS and VLB I techniques. 

An overview of the elements considered in pursuing this analysis Is given 
in Figure 6-9. 

Tracking Configurations (4.1) - In terms of TDAS tracking, the various 
constellations reduce to two situations: frontside and backside satellites. 

Figure 6-10 illustrates configuration options for acquiring TDAS tracking 
data with each technique. 

With BRTS tracking, two-way range and range-rate. (R,R) data are acquired 
from transmissions originating at White Sands (WHS). These are returned 
by bilateration ranging transponders (BRTs) operating at pairs of automated 
ground stations (one pair per TDAS in the minimal configuration). For the 
backside TDAS, transmissions are also relayed via the crosslink which 
impacts tracking accuracy since measurements are affected by uncertainties 
in the frontside satellite orbit. 

With VLBI tracking, the fundamental measurement data type is the difference 
in range (aR) between a signal source and two receivers displaced along a 
known baseline. For the TDAS application, signals received at each auto- 
mated station would be returned to WHS for processing. Typically, one pro- 
cessed pair of VLBI observations per hour per station set (three per TDAS) 
would be needed. Crosslink uncertainties are not a factor due to the 
inherent differencing of common path components. 

Error Modelling (4.2.1) - TDAS OD was assumed to employ batch processing 
of tracking observations taken over a given tracking interval. Orbit pre- 
diction errors beyond the end of tracking were evaluated via a covariance 
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analysis program [16] given a nominal orbit, tracking station locations, 

. tracking schedule and appropriate models of significant tracking error 
sources (measurement noise and bias, station survey, solar pressure and 
tropospheric uncertainties). 

BRTS measurement errors were assumed consistent with TORS OD analyses [18]. 
VLSI measurement errors were considered at two levels: a conservative 
baseline model consistent with values used In studies of a Deep Space Net- 
work (DSN) tracking application to TORS [19] and a more optimistic reduced 
model which reflects achieved capabilities with VLBI stations observing 
celestial radio sources [20]. Station errors were taken to be the same as 
BRTS for the baseline model and an order of magnitude better for the 
reduced model . 

TDAS Tracking Accuracy (4.2.2) - For evaluating performance, the first 
criterion considered was the peak error over a 24 hour continuous predic- 
tion Interval from the end of tracking. Analysis results Indicate, how- 
ever, that since solar pressure is so dominant, OD accuracy can be signi- 
ficantly affected by the choice of reference epoch and duration of the 
tracking Interval and prediction interval. Thus, in both BRTS and VLBI 
analyses efforts were made to Identify multiple combinations (epoch/tracking 
Interval/prediction interval) which, taken together, yield better accuracy 
over a 24 hour period than any one combination alone. 

Table 6-4 gives results for both frontside and backside TDAS OD accuracy 
via BRTS tracking using one 24 hour prediction interval and two concatenated 
12 hour intervals. Multiple segments are effective for the frontside TDAS 
but not the backside which is >2.5:1 higher. The prediction error is due 
primarily to the Impact of the frontside TDAS position uncertainty assumed 
to be 100 M. 

Table 6-5 and Figure 6-11 present the results for TDAS OD accuracy with 
VLBI tracking for both baseline and reduced error models. The reduced 
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TABLE 6-4 

MAXIMUM IDAS POSITION ERROR IN 24 HOUR PREDICTION INTERVAL WITH BRTS TRACKING 


TRACKING 

INTERVAL 

(HOURS) 


150 (SP) 
80 (b) 
40 (s) 

FRONTSIDE 

100*W 




MAXIMUM 
POSITION ERROR 
(M) 

PRIMARY ERROR 
CONTRACTORS* 




BACKSIDE 

98*E 




230 (s') 
175 (sp) 
25 (b) 


220 (s') 
165 (sp) 
25 (b) 


’ LEGEND: SP - SOLAR PRESSURE: B - BRT BIAS; S - STATION SURVEY: 

S*’ - STATION SURVEY f T100 LOCATION (TDAS IOO*V). 

** CONCATENATED INTERVALS BASED ON TRACKING INTERVAL EPOCHS SPACED "12 HOURS APART 
5=1 STANFORD 

I TELECOMMUNICATIONS INC. 



TABLE 6-5 

MAXIMUM TDAS POSITION ERROR IN 24 HOUR PREDICTION INTERVAL WITH VLBI TRACKING 


PREDICTION 

INTERVAL 


VLBI 

TRACKING 

ERROR 

INTERVALS 

HOOEL 

(HOURS) 

BASELINE 

18/36 

BASELINE 
(WITH BIAS 
ESTIMATED) 

18/12 

REDUCED 

6 

REDUCED 


(W. 51 TROPO) 




* TDAS 9 100‘W; STNS: VHS-AGO, WHS-HAW; DATA RATE ■ 1/HR 
** LEGEND: SP - SOUR PRESSURE 
S - STATION SURVEY 
YB - VLBI HEAS. DIAS 
VN - VLBI HEAS. NOISE 
TR - TROPOSHERE 

■I STANFORD 

TELECOMMUNICATIONS INC. 
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model results Illustrate dramatically the benefit of more accurate measure- 
ments since shorter tracking and prediction Intervals can be used to off- 
set the pervasive Influence of solar pressure errors. Additional Improve- 
ment to the 5-8 M level would be possible If better tropospheric correction 
(,< S%) can be achieved. 

Results of the TDAS tracking analysis lead to the following observations: 

• Orbit determination error for backside TDAS satellite tracking 
via BRTS Is significantly higher 2-3:1) than for frontside 
tracking based on two BRT sites per satellite. 

• VLBI tracking offers the potential for significantly Improved 
TDAS orbit accuracy (< 10m) for both front and backside satel- 
lites 

- This method Is not sensitive to cross-link uncertainties, 
but tracking data from three stations would be required 
for each satellite 

- Realization of potential VLBI tracking accuracy requires: 

— Calibration of tracking biases and station location/ 
baseline uncertainties (e.g., using aVLBI data) 

— Judicious selection of epoch (time of day), tracking 
Interval, and prediction Interval combinations. 

• Preliminary results for VLBI tracking with Conus-based sites also 
indicate Increased accuracy, although the shorter baselines may 
reduce the potential improvement compared to that with intercon- 
tinental baselines. 
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6.4 


CONCLUSIONS 


Analysis of the study results leads to the following key findings*. 

• TDAS beacon tracking (FLBT) will satisfy all users In the TDAS 
mission model with position accuracy requirements down to 10 M. 

• Scheduled tracking alternatives (FLST, RLST) can also meet the 
accuracy requirements except at low altitudes where performance 
Is highly sensitive to: 

- Drag uncertainty 

- Frequency of tracking passes, and/or 

- Frequency of navigation data uploads (RLST/only). 

• A two or three satellite TDAS constellation Impacts performance 
as follows: 

- Selecting 2 satellites leads to a tradeoff between coverage 
and accuracy. Increased satellite spacing Improves coverage, 
but a point is reached where performance in high inclination 
orbits begins to degrade (130° spacing appears better than 
162°). 

- Selecting 3 satellites provides full coverage and up to a 
2:1 advantage In navigation accuracy over two satellites. 

• Projected TDAS tracking accuracy requirements (25 M-POS. & 2.5 
MM/SEC-VEL.) can be met with VLB I tracking but not with a mini- 
mal* BRTS configuration. 

6.5 • RECOMMENDATION 


w 


Based on the study results beacon tracking is recommended as the prime 
approach for routine navigation support. Scheduled tracking alterna- 
tives, one-way and two-way, should also be considered for supporting 
user navigation functions, as proposed in Table 6-6. 


2 BRTS sites/TDAS. 
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APPENDIX A 


TDAS MISSION MODEL 


In Task 1 of the TDAS study [lb], a screened baseline of NASA plans was 
used to generate scenarios of experlments/missions for the 1990 - 2005 
time period. An estimated flight schedule was established by first assign- 
ing planned experiments, then candidate experiments and finally opportunity 
experiments. Where planning data for the 1990s was unavailable, generic 
experiments/missions were developed based on trends established In the 
1980* s planning data. Table A-l lists the adopted TDAS mission model and 
corresponding orbit data.* 

Potential navigation accuracy requirements for the various experiments/mis- 
sions were developed from user community survey data, subsequent conversa- 
tions with designated points-of-contar.t and/or Independent estimates. Table 
A-l lists the position and time accuracy requirement;.' developed In each case. 
Except for the TOPEX mission, the most stringent position accuracy require- 
ments Is 10M. Time accuracy requirements were found to be no less than 
1 usee. 

For purposes of comparing user navigation requirements with estimated navi- 
gation performance results, user orbits were divided into six categories 
(A-F). Tables A-2 and A- 3 show the distribution by category and number of 
missions with a given requirement. Some missions are counted more than 
once, since the orbit range spans more than one category (see Table A-l). 


* More detailed information is available in [lb]. 
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TABLE A“1 

TEAS MISSION MODEL (1990 




2005) 



EXPERJMEMT/MISSION 

ALTITUDE 

(to>) 

ORBIT* 

INCLINATION 

(d*9> 

CATEGORY 

ACCURACY 

POSITION 

(») 

REQUIREMENT 

TIME 

(uitc) 


SPACE TELESCOPE 

600 

2B.S 

8 

> 200 

I 


ADVANCED X-RAY ASTROPHYSICS FACILITY 
(AXAF) 

490 

28,5 

6 

> 1000 

J 1 


ASTROPHYSICS GENERIC (AG) • 1.2, 3, 4 

LEO 

28.5 

e 

M 

TBD 


VERY LONG EASELINE RA0I0 
INTERFEROMETER (VLBRI) 

400-5000 

45 

e.c 


1** 

ASTROPHYSICS 

X-RAY OBSERVATORY 

MO 

28.5 

A 

> 200 

>100 

C0SH10 OPTICAL SYSTEM OF MODULAR IMAGING 
C0LLECT0RS/J00 METER THINNED APERATURE 
TELESCOPE 

500 

28 

B 

II 

1 


LARGE OPTICAL/ULTRA VIOLET INTERFEROMETER 

450 

28.5 

B 

* 



ORBITING SUBMILLIMETER TELESCOPE 

1000 

SUN SYNC 

F 

* 

TBD 


LARGE AMBIENT DEPLOYABLE IR TELESCOPE 
(LADIR) 

400-700 

28- SO 

U 

II 

i 100 

SOLAR 

TERRESTRIAL 

SOUR CYCLE ANO DYNAMICS MISSION (SCAOH) 

575 

28 or 98 

S or £ 

> 1000 

1 

SOUR TERRESTRIAL GENERIC (50) • 2 

575 

28 

B 


a 


SOUR TERRESTRIAL GENERIC (SO) - 1,3 

400 

57 

B 

H 



T0PQ5MRHY EXPERIMENT (TOPEX) 

1334 

53.4 

C 

3(0,1 ALT) 

2 ICO 

GLOBAL 

ENVIRONMENTAL GENERIC (EG) « 2,5.6 

1334 

63 

C 

« 

ll 

ENVIRONMENT 

OCEAN SYNTHETIC APERATURE RADAR (OSAR) 

790 

98 

c 

50 

1 


ENVIRONMENTAL GENERIC (EG) • 1,3,4 

790 

98 

£ 

* 

M 


SOIL AND SNOH MOISTURE RESEARCH ANO 
ASSESSMENT MISSION 

400-700 

60-98 

B or E 

100 

w 

RESOURCE 

OBSERVATION 

OPERATIONAL EARTH RESOURCES SATELLITE/ 
ADVANCED UNO OBSERVING SYSTEM 
(OERS/ALOS) 

705 

99 

E 

10 

H 

RESOURCES GENERIC (RG) - 1,2,3 

700 

98 

E 

* 

M 


ADVANCED THERMAL MAPPING APPLICATIONS 
SATELLITE 

620 

97,8 

£ 

SO 

■ 


HAGNETIC FIELD SURVEY B 

300 

97 

0 

30 

« 

METEOROLOGY 

OPERATIONAL METEOROLOGY SATELLITE (2) 

330 

98.7 

F 

500 

_> 100 


SHUTTLE • 1 

185-1110 

28.5-57 

A-C ( 

100 

l (FUTURE) 


SHUTTLE - 2 

135-1110 

70-104 

0-F 

* 

H 

SPACE 

TELEOPERATOR MANEUVERING SYSTEM <TH$) 

tas-iooo 

VARIOUS 

A-F 

N 

TBD 

TRANSPORTATION 

ORBITAL TRANSFER VEHICLE (OTV) 

LEO-GEO' 

VARIOUS 

A-F 

M 

M 


HEAVY LIFT UUNCH VEHICLE (KLLV) 

200-500 

VARIOUS 

A,B or O.E 

M 

w 


MANNED ORBITAL TRANSFER VEHICLE (HOTV) 

LEO-GEO 

VARIOUS 

A-F 

M 

H 


SPACE STATION/PUTFORM (2) 

400 

28. S 

s 

10-30 

1 


* LEO • Low Altitude Earth Orbit 
GEO - Geostationary Earth Orbit 

** For Clock Synchronization Only. Needs Hydrogen MASER Stability, 
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TABLE A-3 

USER TIKE ACCURACY REQUIREMENTS AND ORBIT DISTRIBUTION 

IN TDAS MISSION MODEL 
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APPENDIX 3 


TRACKING DATA MEASUREMENT ERRORS: RANDOM COMPONENT 

Random errors in the range and Doppler measurement data arise due to noise 
in the tracking link. This appendix estimates these errors for the three 
one-way tracking alternatives based on the signal definitions and assump- 
tions in Section 3. Parametric results are presented where applicable to 
Indicate the sensitivity to design parameters; e.g., TDAS EIRP and user G/T 
or user EIRP. The results are presented as follows: 

B.l Thermal Noise in Range Measurements 
8.2 Thermal Noise in Range-Rate Measurements 
B.3 Phase Noise in Range-Rate Measurements 

B.l THERMAL NOISE IN RANGE MEASUREMENTS 


The measurement of range in the systems considered here relies upon the 
acquisition and tracking of a pseudonoise (PN) ranging code. In all the 
one-way navigation modes, the receivers exploit the correlation properties 
of PN codes to produce an error signal which directs the timing correction. 
The principal figure of merit is the variance of the timing error as a 
function of the input ratio of signal power to noise density, C/N 0 . 


B.1,1 Forward Link Tracking 


The Motorola 2nd generation TDRSS user transponder design [14] is assumed 
to apply to the users to be supported in the TDAS era. The design's code 
tracking loop is a digital Implementation of a 2nd order type II tau-dither 
loop; a rate aided tracking loop configuration is also present, providing 
a 1st order tau-dither loop with frequency tracking from the carrier loop. 
Motorola expresses the variance of the tracking jitter due to thermal 
noise as 



(chips) 2 , 
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where Bl is the closed loop bandwidth of the loop, either 1 Hz for the 
independent tracking mode or .125 Hz with rate aiding, and Bjf is the IF 
filter bandwidth of 4 kHz. This result can be seen to stem from Hartmann's 
analysis [25] of the tau-dither loop where the incoming PN coded signal is 
alternately correlated with early and late versions of a locally generated 
code replica. 

Based on the forward link budget in Figure 3-1, Figure B-l shows the RMS 
range error in meters as a function of the user G/T, parameterized by the 
TDAS EIRP. The mean square range error cr^ is derived from the mean square 
timing error a ^ by the simple conversion 

a (c T c a T ) 2 = (d c o T ) 2 (meters) 2 , 

where c is the speed of light and T c is, again, the chip duration, assumed 
here to be 1/3.0778 Mcps. From these curves the forward link range errors 
due to thermal noise are extracted for a user with G/T = -27 dB/°K for 
both FLBT and FLST modes. 

8.1.2 Return Link Tracking 


For return link tracking, a ground receiver incorporating the Harris wide 
dynamics demodulator (WDD) is envisioned. The WDD uses a delay-locked loop 
to track the PN coded signal: early and later versions of the local PN 

code are correlated simultaneously with the input PN code sequence. An 
analysis of the exact implementation used in the WDD is not readily availa- 
ble, but work by Simon [26] provides results applicable to a delay-locked 
loop where the bandpass arm filters have bandwidths on the order of the 
data rate. A linear analysis for the mean-square timing jitter, assuming 
a large equivalent loop signal-to-noise ratio yields 


a 2 s „ik_ 

°T 2 C/N 0 


jj- (chips) 2 . 
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where Bl is the single-sided tracking loop bandwidth and Sl is the squaring 
loss of the delay-locked loop, Simon shows that 


SL 




Kd + k l 


ZBarn/Rs 
. R d D m . 


where 


R s = data symbol rate = 1/T S , 

B arm = 2-sided arm filter noise bandwidth 
R d 3 symbol SNR = jp T s , 




ho 


f s d (f))H(f)|2 df, 


Kd 


ha 


S s d(f ) |H(f )| 4 df 


ho 


/ S d (f))H(f)|2 df 


XL 




/ |H(f) | 4 df 

“CD 

7 |H(f)|2 df 


S d (f) a baseband power spectral density of the data, 
and H(f) = transfer function of the arm filter. 


In keeping with the WDD parameters, it is assumed that B^ = 4 Hz in the 
tracking mode and that 1-pole arm filters are used. The symbol rate used 
for the calculation is assumed to be the maximum achievable data rate for 
a given user spacecraft El RP , as detailed in Section 3.2.2 and shown in 
Figure 3-5. With the anticipated rate 1/2 convolutional coding, the 
symbol rate R s is found as twice the achievable data rate. The tracking 
filters used in the WDD are chosen on the basis of symbol rate to optimally 
band! imit the input noise. The 3 dB bandwidth f c for the WDD is thus 
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8 kHz for all R s < 4 ksps 

39 kHz for all R s < 25 ksps 

300 kHz for all R s < 300 ksps 


where Sip « tt f c for a 1-pole filter, H(f) 
for a 1-pole filter is found as: 


3 1 
1+jf/fc ' 


kl 



1 

1 


The parameter K|_ 


The other parameters in Simon's expression depend on the data modulation 
format as well as the arm filter transfer function. Assuming equiprobable 
independent transmitted symbols — an assumption not strictly valid with 
coding — the power density for the NRZ signalling format is given by 


S d (f) 3 T s 


sin 2 (rfT s ) 

(irfT s )2 


y 


for biphase (Manchester) coding, ( f ) 
that 


sin4(mfT s /2 ) 

S (irfT s /2) * 


It can be shown 


°m 


1 L- ( i« e -2irR * ) for NRZ coding 

2irR ' 

1 - ~-r (S-Ae" 71 ^ + e" 2lT R ) for biphase coding 


KD 



4irR* 


for NRZ coding 


>4(3-hrR‘ + (3+27rR‘)e-fr r R'] 

4ttR ‘ 


for biphase coding 


where R’ = f c T s . 
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Figure B-2 shows the above results expressed as an RMS range error in 
meters versus the user EIRP corresponding to the return link budget in 
Figure 3-5. The discontinuities appearing on the graph result from the 
change In arm filter bandwidth with data rate. A user EIRP of 2 dBW is 
assumed as a baseline value. 

B.2 THERMAL NOISE IN RANGE-RATE MEASUREMENTS 

Range-rate is estimated on the basis of a signal's observed doppler shift 
in frequency due to the relative motion between the transmitter and 
receiver. A carrier-tracking loop in the receiver provides the means to 
extract the phase <t> (in radians) of the received signal and thus the 
received frequency, since <p 3 2 w ft. By counting positive-going zero- 
crossings, N(t), over the averaging time T av , an estimate r of the range 
rate R is derived as • 

i, » x [ N ( t+T av) " N(t) = r + _L Mt+Tay) ~ Mt) f 
L iav -1 2n T7v ~ 

where 6<j>(t) represents noise in the phase measurement and \ is the carrier 
wavelength, assumed to be known a priori for one-way navigation. Assuming 
T av is long enough so that the phase samples are independent and identi- 
cally distributed, the RMS range-rate error is given by 

= y oi *<t+T aV ) + <(1) - ^ »♦ 

Here is the RMS phase error in radians and one may identify as the 
RMS doppler phase error. The performance of the carrier-tracking loop in 
the receiver thus determines the range-rate error due to thermal noise. 

In the results presented here, the doppler averaging time T av is assumed 
to be 1 second. 
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B.2.1 Forward Link Tracking 


The Motorola 2nd generation TDRSS user transponder's carrier-tracking loop 
Is a 2nd order type II modified Costas loop. For Costas loops In general, 
the mean square phase error Is given by [27, 28]: 

»♦ ■ (rad,iins)2 • 

where Bl Is the loop noise bandwidth and Sl Is the Costas loop squaring 
loss. The phase of a Costas loop, however, Is twice the received signal 
phase, so the thermal noise phase jitter 0 $ that Is of concern here Is 

2 - 8 L 1 
% " C/N 0 * Sl 

For the Motorola design, Bl a 40 Hz and $l = .73. With these parameters, 
Figure B-3 shows the RMS one-way range-rate error anticipated for forward 
link tracking assuming the link budget of Figure 3-1 and the various possi- 
ble TDAS EIRP values. As before, a user G/T value of -27 dB/°K Is used as 
a baseline. 


3.2.2 Return Link Tracking 


The wide dynamics demodulator's 1st order carrier-tracking loop, shown 
In Figure B-4a, has been analyzed by Weinberg [29], yielding the following 
expression for the received mean square phase error: 


a 


2 

<P 


„ 8Us 
a B E|j/N 0 


1 

1 " 2 b L t s 
8 J 


( radians) 2 


Here 

Ts 

ik 

No 


symbol duration, 

symbol SNR 3 T s , 

N o 

i - 2 « e ' Eb/N ° ■ 


P 

f ; 


n > 

* « 
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FIGURE B-4: RLST: CARRIER TRACKING LOOP' ERROR DUE TO THERMAL NOISE 


NmowwK) 

AVIM01N4 

nun 



mmiUMwa 


PERFORMANCE MEASURES 

RESULTS* 

WIDE DYNAMICS DEMODULATOR (PLL) 
PHASE ERROR, - y^T 

109 MRAD (6.2*) 

DOPPLER PHASE ERROR - vT 

154 MRAD. (8.8*) 

RANGE-RATE ERROR - V? o, 

(1-WAY) 2,7 Tav 

3.5 MM/SEC 


- BER - 10-5, 3 DO MARGIN, 5 DB CODING GAIN 

- Eq/Ho - 4.3 DB (WORST CASE) - (C/N 0 )/SYMBOL RATE 

K - .032, A FUNCTION OF E B /N 0 , LOOP BANDWIDTH AND SYMBOL RATE 

- T av - 1 SEC, DOPPLER AVERAGING TIME 

A - SMA WAVELENGTH (MM) 
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b - l - Q( n £b/N 0 “)i 

and Q{x) » complementary Gaussian distribution function 
■ Jtn / e" u2 / 2 du 

In the WDD, the loop bandwidth Is adjusted proportional to the symbol rate; 
for this calculation, It Is assumed that B[Js % <03. If worst-case opera- 
tion Is assumed and the maximum achievable data rate Is used for all values 
of user EIRP, then Figure 3-5 shows that Eb/N 0 Is constant at 4.3 dB-Hz. 
Consequently, the carrier tracking error Is constant, as summarized In 
Figure B-4b. 

B.3 PHASE NOISE IN RANGE-RATE MEASUREMENTS 

Additive thermal noise Introduced by the channel Is not the sole source 
of random fluctuations that affect range-rate measurement precision. A 
system's transmitting and receiving oscillators, system mixing chains, AM 
to PM conversion In the channel, spurious vibrations, etc. all Introduce 
phase noise components Into the signal which degrade the carrier-tracking 
loop performance of the receivers. To estimate the Impact of phase noise 
on one-way navigation, consider first the doppler measurement process: 
the phase value <p extracted by the carrier-tracking loop Is sampled at 

time t-T av and at time t and then differenced. That difference A<t> Is 

• • 

related to the range-rate R as R 3 ( x/2irT av ) a</> . Since a* ° <t(t) - 
<p( t - T av ), the transfer function describing the operation Is simply 
F(f) 3 l-e"J27rfT aVt if the Input to the doppler extractor has single-sided 
phase noise spectral density S^(f), then the spectral density of a* is 
given by 

W f > = |F(f) | 2 S 0 (f) 3 4 sin 2 (irfT av ) S*(f). 

Since o 2 ^ => 2 f (f) df, the input spectral phase noise density S<$(f) 


VI-8—11 


contributed by the system can be seen to determine the range-rate mean 
square error due to phase noise. 

In this analysis, only the effects of phase noise In the system oscillators 
and VCO's as propagated through the system mixing chain are considered, 

The key to the analysis, as discussed In [30], lies In characterizing a 
carrier-tracking loop by Its closed loop transfer function H(f). A linear 
analysis of the carrier loop shows that the loop's VCO output tracks the 
carrier with noise contributions from two sources. One Is the Input phase 
jitter lowpass filtered by H(f), the other, the VCO's own phase jitter 
hlghpass filtered by l-H(f). In other words, If the Input phase noise 
spectral density Is S 0 (f) and the phase noise spectral density of the VCO 
itself Is S 0 ' ( f ) » then the coherent reference provided by the tracking 
loop has phase noise spectral density S 0 (f) « | H( f ) | 2 S$(f) + |l-H(f)| 2 Sg -(f). 
The simple analysis explained here has been validated by extending It to 
model the two-way doppler tracking of TDRSS ; the results are consistent 
with other estimates of the doppler measurement error due to phase noise. 

The forward and return links are modeled In terms of the ground and user 
oscillators, the tracking loops and their VCO's, and the frequency multi- 
plies comprising the mixing chains. It Is assumed that TDAS satellites 
will operate In the same fashion as do TDRSS satellites; a ground-generated 
pilot tone is tracked on-board the TDAS satellite to provide the reference 
tone for frequency-transl tlon of the relayed signals, TDRSS uses a K u -band 
space-to-ground link; In TDAS, both K u - and K a -band frequencies are being 
considered. Figure B-5 shows the models used here for the forward and 
return links, respectively. Both the user and ground station's oscillators 
are assumed to be 5 MHz standards; the frequency multiplies required to 
support return and forward link service and a K u - or K a -band space-to-ground 
link are shown in the diagrams. The Ku-band uplink, downlink, and pilot 
frequencies are taken from TDRSS; those for K a -band are chosen arbitrarily. 
S-band forward link service centered at 2106.406 MHz (as In TDRSS) and an 
assumed IF frequency at the user spacecraft of 9.5 MHz yield a fixed 
multiply value at the user oscillator, as shown In Figure B-5a. S-band 
return link service centered at 2287.5 MHz similarly requires a fixed 
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FIGURE B-5: PHASE NOISE MODELS FOR TDAS FORWARD AND RETURN LINKS 
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multiply value at the user oscillator In Figure B~5b. For return link 
tracking, an IF of 240 kHz Into the ground carrier-tracking loop is assumed. 

The generation of high-frequency signals from the ground and user frequency 
standards Is pessimistically assumed to require frequency synthesizers and 
hence tracking loops and their VCO's, Noise from VCO's, of course, is 
also introduced at the pilot tracking loop on-board the TDAS satellite 
and at the ground and user receivers' tracking loops. For simplicity, 
all the tracking loops are assumed to be 2nd order with closed loop trans- 
fer function 

s w n 2 +j/Tw n w 
{ 1 w^V jyFw n w - w* 


where w 3 2rrf, w n Is the loop natural frequency, and the single-sided loop 

noise bandwidth Bi % P... . Each of the tracking loops may thus be para- 

4 iz 

meterlzed by its loop bandwidth. The following values are chosen on the 
basis of present TDRSS practice, the Motorola 2nd generation TDRSS user 
transponder design, and arbitrary selection: 


Hi ( f ) (ground frequency system): 
HgCf) (TDAS frequency system): 
H3(f) (user frequency system): 
Hq(f) (user tracking loop): 

Hg(f) (ground tracking loop): 


Bl 5 700 Hz, 
B L = 262 Hz, 
B L 3 700 Hz, 
Bl = 40 Hz, 
B L 3 70 Hz. 


To complete the model, the single-sided phase noise spectral densities of 
the ground and user frequency standards and the VCO's are approximated 
as shown in Figure B-6. As in TDRSS, the ground frequency standard is 
assumed to be an HP 5061A+004 cesium oscillator with its phase noise spectral 
density derived from the manufacturer's specifications [31]. The user 
frequency standard is assumed to be a precision quartz oscillator; its 
phase noise spectrum is estimated from the specifications for the FTS 
1150 quartz standard [32]. Note that Ulile the cesium oscillator provides 
superior long term stability, its short term phase fluctuations are more 
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severe than those of the quartz oscillator. In all the figures, the 
spectral density is assumed to extend for offset frequencies beyond those 
shown in the same manner as immediately before the end of the depicted 
range. Only for the VCO spectral density will these high offset frequencies 
be significant since only the VCO output Is highpass filtered by transfer 
functions of the form l-H(f). 

8.3.1 Forward Link Tracking 

Figure B~7a explicitly shows the various filtering operations to which the 
phase processes <f>cESIUM> VCO » and tfxTAL in the forward link are subjected. 
The ground-to-TDAS delay n and TDAS-to-user delay t 2 are accounted 
for by the appropriate phase shifts; n and t 2 are approximated here as .14 
seconds, corresponding to distances of 42000 kilometers, Inspection 
of the diagram allows the transfer functions associated with the ground, 
TDAS, arifl user segments to be identified separately. This yields the 
component single-sided phase noise spectral densities at the input to the 
user's doppler extractor; 

W f) = S *CESIUM (f)|H4<f)Hl(f)e ' jW<tl + T 2 , (K 2 -K3H 2 {f))|2 

- S* vco (f)|(l-H 1 (f))H 4 (f)e-.MT 1 + T 2 )(K 2 .K 3 H 2 (f))|2 

' [VsiUH (f,|Hl(f) l 2 + \cO (f)|1 - Hl(f)|2 ] X 

|H 4 (f)|2[ld;- 2K 2 K3Re(H 2 (f)> + k|| H 2 ( f ) J z] ; 

s W f) = \co |(1 - H2(f,)K3 H4(f)|2; 

S «USEP. (f) = [ S *XTAL (f)|H3(f)|2 + 3 * vc0 (O|l-H3(f)| 2 ](419.38)2|H 4 (f) 

+ s W f)|l ' H4<f)|2i 

where Re{ • > denotes "the real part of {•}", 
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FIGURE B-7s DERIVATION OF TRANSMISSION LINK TRANSFER FUNCTIONS 
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The phase noise variances corresponding to each of these densities Is 
obtained as simply a ^x 3 2 I S $x (f) df ; the range-rate variance due to 
each component is similarly obtained as 

C R 2 3 2 I s *x ^ 4 si ° 2 ( 7 rfT av) d ^ 

The results obtained from numerical integration are reported in Table 3-3. 


3.3.2 


RETURN LINK TRACKING 


The necessary transfer functions are found from Figure B-7b in the same 
way as in the preceding section. As previously, the single-sided phase 
noise spectral density at the input to the doppler extractor may be 
broken into its component elements according to the phase noise sources: 


S 

S 

S 


*USER^ 

*TDAS (f) 

*GND (f) 


S *XTAL (f) ^ 3 (f) l 2 + S W f) 11 _H 3 (f)|2 ( 457 ' 5 ) 2 ll H s (f) S 2 ’’ 

» « 


s +vco (f)i(1 ' H2<t),l<4 H5(f)|2; 
S *CEsW H l (f) l 2 + Vo |1 ' Hl(f)| 




l< 4 |H 2 (f )| 2 + K 5 2 - K 5 K 4 C 2 Re(H 2 (f)}cos 2 wtl 
+ 2 Im{H 2 (f)}sin 2 wt 1 ] 


+ S % 0 (f)|l-H 5 (f')| 2 i 


where Im {•} denotes "the imaginary part of {•}". 


Again, the range-rate variance a« due to each component is 

R x 

00 

= 2 / 4 sin 2 {irfT av )df. 

I 'X • J Q 


The results from numerical integration are shown in Table 3-5. 
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APPENDIX C 


TDAS TRACKING ANALYSIS - SUPPLEMENTARY RESULTS 
C.l INTRODUCTION 

An evaluation of the potential tracking accuracy for TDAS satellites was 
made based on BRTS and VLB I tracking techniques. Section 4 presented a 
brief description of each technique, possible tracking configurations to 
support frontside and backside satellites, and the error analysis modeling 
approach and major results. This appendix provides supplementary results 
for BRTS tracking in Section C.2 and VLSI tracking In Section C.3. Relevant 
modelling data used In the analysis are given In Section C = 4, 

C.2 BRTS RESULTS 

TDAS position errors in the prediction interval are presented In Figures 
C-l through C-4 as a function of BRTS tracking epoch and tracking Interval. 

Figure C-L (Frontside TDAS) illustrates the sensitivity of TDAS position 
errors In the prediction Interval as a function of epoch time. In all cases, 
the minimum error bound is about 100 m, regardless of epoch time. However, 
the maximum error bound Is cyclic with roughly a 12 hour period, and the 
peak to peak variation diminishes with longer tracking Intervals. To mini- 
mize the position error In the 12 hour prediction Interval, an epoch Is 
chosen at the valley of the maximum error bound. Note that the valley 
epochs shirt as the tracking Interval changes. Of the four tracking inter- 
vals shown (18, 24, 30, and 36 hours), the 18 hour case has the lowest val- 
leys (110 m), which occur near 18 H 30 M GMT and S H 30 M GMT. The 24, 30, and 
36 hour cases all have valleys' on the maximum error bound greater than 125 m, 
which occur at epochs different from the 18 hour case. For example, the val- 
leys for the 24 hour case occur near 0 H 30 M anc j 12 H 30 M gmt, w ith a magnitude 
of 140 m. For the 24 hour prediction Interval, the valleys In the maximum 
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. FIGURE 04 

TDAS MAXIMUM POSITION ERROR CONTRIBUTORS IN 12 HOUR PREDICTION 
INTERVAL VS. BRTS TRACKING INTERVAL 
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error bound are all above 140 m for all tracking Intervals, Thus, for the 
frontside TDAS (at 100°W), the best choice Is 18 hour tracking, 12 hour pre- 
diction, and epoch of 18 H 30 M , resulting In TDAS la position errors between 
100 and 110 m. 

Figure C-2 demonstrates the same relationships for the backside TDAS at 
98°E, for which the best choice Is 36 hour tracking, 12 hour prediction, 
and epoch of 12^30^1, resulting In TDAS la position errors between 210 and 
250 m. However, note that the maximum error bound Is relatively flat for 
the 36 hour tracking Interval and remains below 300 m over all epochs. 

Thus, the epoch Is not as critical for the 36 hour tracking case as It Is 
for 18, 24, and 30 hour tracking cases In which the maximum errors fluctuate 
widely (exceeding 400 m). 

Figure C-3 Illustrates the TDAS position errors from the end of tracking to 
the end of prediction for a fixed epoch time. Figure C-3a Is for TDAS at 
1Q0°W, with a tracking configuration using transponders at White Sands (WHS) 
and Santiago (AGO). WHS-WHS denotes a roundtrip path from WHS transmitter- 
to TDAS to WHS ground transponder and return to TDAS to WHS receiver. Simi- 
larly, WHS-AGO denotes a roundtrip path from WHS transmitter to TDAS to AGO 
transponder and return to TDAS to WHS receiver. The TDAS position errors 
are plotted for 18, 24, 30, and 36 hour tracking intervals, all for 24 hour 
prediction Intervals, based on an epoch time of 18 H 30 M GMT. Again, note the 
best prediction Interval occurs during the 12 hours after the end of the 18 
hour tracking. Figure C-3b is for TDAS at 98°E with transponders on the 
100°W TDAS (T100) and at Diego Garcia (DGA). T100-T100 represents a round- 
trip path from T100 to TDAS (98°E) and return to T100.* T100-0GA represents 
a roundtrip path from T100 to TDAS to DGA transponder and return via TDAS 
to T100.* The best prediction interval for the backside TDAS is during the 
12 hours after the end of the 36 hour tracking. 


For error analysis with existing ORAN capabilities the WHS-T100 link 
was not modelled. Instead, separate WHS-T100 path delay measurements 
were assumed to be subtracted from backside two-way BRTS measurements 
made at WHS. ORAN was set up to use T100 as the terminal location 
(analogous to WHS) with station errors equal to the position accuracy 
for frontside TDAS tracking. 
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Figure C-4 presents the TDAS maximum error contributors in the 12 hour pre- 
diction Intervals for 18, 24, 30, and 36 hour tracking intervals. Note the 
position error scale is a log scale. Figure C-4A is for TDAS at 1Q0°W, and 
shows solar pressure, VLBI measurement bias, and station survey errors as 
dominant, with magnitudes greater than 40 m. Figure C-4B is for TDAS at 
98°E, and shows station errors and solar pressure uncertainties as the 
domin. . ‘ r >rror sources, with each contributing well above 100 m. Note that 
the station error reflects uncertainty due to the frontside TDAS (T100) 
tracking error assumed to be a constant 105 m (H a 26 m, C ■ 28 m, L * 98 mj . 


C.3 


VLBI TRACKING 


This section presents TDAS tracking accuracy results based on the VLBI 
tracking technique. Subsections C.3.1 and C.3. 2 discuss results generated 
for the baseline and reduced error models defined in Section 4 (Table 4-1). 
Subsection C.3. 3 presents some preliminary results for frontside TDAS 
tracking using VLBI stations located only In CONUS. These three subsec- 
tions show the effects of tracking Interval, epoch time, prediction 
Interval, VLBI bias estimation, error sources, and data rate on TDAS 
uncertainty. 

C.3.1 Baseline Error Model Results 


A composite plot of TDAS lo position uncertainty during the prediction 
Interval is shown in figure C-5 for several treking Intervals - 3, 6, 12, 
18, 24, 30, 36, and 42 hrs, all using a fixed epoch time of 18^30 M GMT. 
Figure C-5A shows the error profile without VLBI bias estimation. For the 
3, 6, 12, 24, 30, and 42 hr tracking, the position uncertainties rise 
above 100 m after a few hours in the prediction interval. Only for the 
18 and 36 hr tracking does the position uncertainty remain below 100 m 
for several hours In the prediction Interval. As will be shown later, 
concatenation of several 18 and 36 hr tracking segments may provide 24 
hours of prediction below the 100 m level. 

Figure C-5B shows the same cases as C-5A except VLBI bias Is estimated. 

The best choices of prediction Intervals are those following the 12 and 
18 hr tracking. These may be concatenated to provide 24 hours of prediction 
below the 100 m level. Note that the position error Is close to 200 m at 
the end of the 36 hr tracking Interval. 
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Figures C-6 and C-7 display the TDAS maximum error contributors for various 
prediction and tracking Intervals with and without VLSI bias estimation. 

As mentioned above, the 18 and 36 hr tracking intervals are the best choices 
for either 6 or 12 hr prediction Intervals. In Figure C-6A for 6 hr tracking 
note that VLBI bias station survey and noise errors are dominant, while for 
long tracking Intervals (24, 30, and 36 hrs) solar pressure error is dominant. 
For Intermediate tracking intervals (12, 18 hrs) VLBI bias, solar pressure, 
and station survey errors are dominant. As the prediction interval increases, 
however, solar pressure ultimately dominates. that in all of these 

cases, the gravity and tropospheric effects are relatively insignificant 
(less than 20 m). 

Figure C-7 is based on VLBI tracking with bias estimation. As r.tfted before, 
the 12 and 18 hr tracking intervals are the best eftciees foF 6 or 

12 hr prediction Intervals. With VLBI bias estimatedi ctfcar related errors 
(station survey and tropospheric errors) are also reduced whereas noise and 
solar pressure error increase, particularly noise. The overall result, 
however, is a decrease in the total error. 

Figure C-8 shows the sensitivity of the TDAS position error profile to epoch 
time of day at the beginning of an 18 hr tracking interval with VLBI bias 
estimation. For each prediction Interval (1, 3, or 6 hr), the minimum and 
maximum error bounds are plotted. The best epoch choices are 6^30^ GMT 
(TDAS local midnight) and 18 H 30 M GMT (TDAS local noon). The minimum error 
bound Is in the 48-55 m range for all three prediction Intervals. The 
maximum error bound Increases dramatically from 1 hr prediction to 6 hr 
prediction when the epoch is not near 6^30^ or 18 ^ 30 ^, For these best 
epochs, there is little change in maximum error bound as the prediction 
interval Increases from 1 to 6 hr. 
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FIGURE C-6 

TDAS MAXIMUM POSITION ERROR CONTRIBUTIONS VS VLSI TRACKING INTERVAL 

(WITHOUT VLB I BIAS ESTIMATION) 
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C.3.1.1 Impact of TDAS Location , Figure C-9 illustrates that TDAS position 
uncertainty is insensitive to TDAS geographic location. Four prediction 
profiles are plotted, based on TDAS longitudes of IOO°W, 98°E, 171°W, and 
4i°W. Both the 18 and 36 hr tracking cases are shown. The small differences 
among each set of four profiles can be attributed to geometrical differences 
In the tracking station configurations. 

C.3.1.2 Effect of VLBI Data Rate . Table C-l is a tabulation of minimum 
and maximum TDAS position uncertainties as a function of data rate, tracking 
Interval, and prediction interval for an epoch of 18 ^ 30 ^ GMT. For a 6 hr 
tracking interval, the effect of data rate is dramatic with the maximum 
error increasing by more than 2:1 as the data period Increases from 60 sec 
to 3600 sec. For a 12 hr tracking Interval, the factor is about 20%. For 
longer tracking intervals (18, 24, and 36 hours), there are virtually no 
differences in position errors for data periods of 60, 1800, or 3600 sec. 

This has significant operational importance since VLBI measurements data 
collection and processing need occur no more frequently than once per hour. 

C.3.1.3 Prediction Interval Concatenation . Figures C-10 and C- II demon- 
strate how it is possible to concatenate several prediction segments to 
achieve better TDAS position performance over a 24 hr period compared with 
a single segment. The 18 hr tracking segments are shown in Figure C-10A. 

Two start times are Indicated - S/C local midnight and S/C local noon on 
3/22. At the end of tracking of these two cases, the 9 hr prediction 
profiles may be used, i.e., on 3/23 the interval between 0^30^ and 9^30^ 
and the second interval between 12 H 30 M and 21 H 30 M . Gaps from 9 H 30 M to 
12 ^ 3 C)M and from 21^30 M to 24^30^ may be filled using prediction intervals 
based on 36 hr tracking starting at S/C local noon and S/C local midnight 
on 3/21 and 3/22, respectively (see Figure C-10B). When all these segments 
are concatenated in Figure C-11A, the result is a 24 hr prediction Interval 
in which the TDAS position error remains below 85 m. 

The same procedure may be used when the VL8I bias is estimated (see Figure 
C-11B). By concatenating multiple prediction segments following 18 and 


VI-C-14 


FIGURE C-9 

HAS POSITION UNCERTAINTY IN PREDICTION j. . ,<VAL VS TDAS LOCATION 


ORIGINAL PAGE ffl 
OF POOR QUALITY 



STANFORD 

TELECOMMUNICATIONS INC. 




ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE C-l 

MIN/MAX POSITION UNCERTAINTY (M) IN PREDICTION INTERVAL 
VS. VLB I DATA RATE AND TRACKING INTERVAL* 


PREDICTION 

TRACKING 

INTERVAL 

INTERVAL 

(HOURS) 

(HOURS) 


6 


12 

6 

18 


24 


36 


„ 

6 


12 

12 

18 


24 


36 


WITH BASELINE ERROR MODEL 


VLBI DATA RATE (SEC) 

60 I 1800 


66/180 

49/110 

53/63 

53/131 

71/92 


66/481 

49/180 

53/131 

53/219 

71/131 




92/437 

52/127 

53/65 

53/128 

71/90 


92/1108 

52/210 

53/130 



TRACKING CONFIGURATION: TDAS 9 1Q0°W; STNS: WHS-AGO, WHS-HAW 

EPOCH: 86/3/21, la^O** 
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FIGURE C-10 

TDAS POSITION UNCERTAINTY VS TIME AND VLBI TRACKING EPOCH 


A) 18 HOUR TRACKING INTERVAL 



B) 38 HOUR TRACKING INTERVAL 
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FIGURE C-ll 

IDAS' POSITION UNCERTAINTY IN 24 HOUR PREDICTION INTERVAL 
(BASED ON MULTIPLE VLB I TRACKING SEGMENTS) 

(A) BASELINE ERROR MODEL 
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12 hr tracking Intervals, an error profile that remains below 70 m can be 
obtained. Also shown In Figures C-11A and B are the radical, cross track, 
and along track components of TDAS position error. In both cases, the 
along track component dominates. 

C.3.2 Reduced Error Model Results 

TDAS position uncertainties computed from the baseline and reduced error 
models are compared In Figure C- 12. Four tracking Intervals were used - 
1, 3, 6, and 12 hrs. For the baseline model, longer tracking intervals 
produce better prediction profiles, as observed In Subsection C.3.I. 

However, for the reduced error model, shorter tracking Intervals appear 
advantageous. A 6 hr tracking Interval followed by a 1 hr prediction 
interval seems to be an optimal choice. As will be demonstrated later, 

1 hr segments may be concatenated to form a 24 hr prediction 'Interval 
with TDAS position errors between 10 and 12 m, a significant Improvement 
over the baseline model. 

TDAS maximum position error contributors are shown In Figure C- 13 for the 
1, 3, 6, and 12 hr tracking Intervals. Note tnat the tropospheric error 
contribution Is significant here, but decreases, while the solar pressure 
contribution Increases with tracking interval. A crossover between these 
two error sources occurs around a 6 hr tracking Interval. All the other 
error contributors are relatively small In this region. 

C.3.2.1 Im pact of Tracking Data Rate . Table C-2 lists minimum and maxi- 
mum TDAS position uncertainties as a function of data rate, tracking 
interval, and prediction Interval. For 1 hr tracking, the 3600 sec period 
(1 measurement per hr) Is Inadequate for orbit estimation (number of measure- 
ments less than the number of unknowns being solved for). For 3 hr tracking 
the 1/hr rate produces a maximum error about twice that of the 1/min rato. 

For 6 or 12 hr tracking, changing the data rate from 1/mln to 1/hr has 
little impact on the results, 
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TABLE C-2 

MIN/MAX POSITION UNCERTAINTY (M) IN PREDICTION INTERVAL 
VS. VLB I DATA RATE AND TRACKING INTERVAL* 


WITH REDUCED ERROR MODEL 


PREDICTION 

INTERVAL 

(HOURS) 

TRACKING 

INTERVAL 

(HOURS) 

VLBI DATA RATE (SEC) 

60 | 1800 | 3600 


1 

96/114 

333/398 

X 

1 

3 

14/16 

19/23 

23/29 


6 

11. S/18.7 

11.6/18.9 

11.7/19.3 


12 

11.6/15.1 

11.4/14.7 

11.2/14/4 


1 

96/212 

333/740 

X 


3 

14/27 

19/43 

. 23/56 

1 *3 

6 

11.5/18.7 

11.6/18.9 

11.7/19.3 


12 

11.6/28 .2 

11.4/27.6 

11.2/27.1 


Vt-C-22 


tr:J 














C.3,2.2 Prediction Interval Concatenation , Figure C-14 illustrates how 
a series of 1 hr prediction segments may be concatenated to form a 24 hr 
prediction profile. Each segment is based on 6 hours of tracking data. 

In (A), ti'4 dominant error source is the tropospheric error (15%),* and 
the profile Is relatively flat over 24 hours with a mean TDAS position 
error near 11.5 m. In (B), the effect of a possible reduction in tropo- 
spheric error to 5% is shown with TDAS position error reduced to the 5-8 m 
level. Solar pressure is the dominant contributor as is evident from the 
12 hr cycle. 

C.3.3 CONUS-Based VLBI Tracking Results 

Some preliminary results on TDAS position uncertainty based on CONUS tracking 
sites rather than global sites are presented in Figures C- IS and C- 16. The 
CONUS sites used were GSFC, HSN (Houston), and SUN (Sunnyvale). Figure C-15A 
shows the prediction error profiles generated with the baseline model for 
6-48 hr tracking Intervals in 6 hr steps. As observed in Section C.3.1, 
longer intervals give better performance than shorter Intervals with measure- 
ment and station survey errors at the baseline model levels. Figure C-15B 
indicates that the maximum TDAS position error during a 6 hr prediction 
interval is dominated by either solar pressure or VLBI bias for the tracking 
Intervals considered. 

Figure C- 16 demonstrates how several prediction Interval segments can be 
concatenated to form a 24 hour prediction error profile. In (A), with 
segments based on 24 and 48 hr tracking intervals, the TDAS position uncer- 
tainty ranges between 80 and 125 m. In (B), VLBI bias Is estimated, and 
the TDAS position uncertainty decreases to the 50-75 m level. 

C. 4 MODELLING DATA 

Tables C-3 and C-4 list various modelling data used in the analysis. 


A 15% uncertainty in tropospheric delay modelling was assumed in the 
VLBI error model (Table 4-1). 
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FIGURE C-14 

IDAS POSITION UNCERTAINTY Iil 21 HOUR PREDICTION INTERVAL 
(BASED ON MULTIPLE VLBI TRACKING SEGMENTS) 
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FIGURE 015 

TDAS POSITION UNCERTAINTY IN PREDICTION INTERVAL AND 
ERROR CONTRIBUTORS FOR CONUS BASED VLBI TRACKING 


(A) POSITION ERROR VS. TRACKING INTERVAL 
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APPENDIX D 


USER NAVIGATION PERFORMANCE RESULTS - SEQUENTIAL DATA PROCESSING 

TDAS user navigation performance was evaluated In terms of potential OD/TD 

♦ 

accuracy for the three one-way tracking alternatives defined In Section 
2. The objective was to assess several user orbit types and the Impact 
of two and three satellite TDAS constellations, the -tracking schedule and 
two algorithms for tracking data processing (sequential and sliding batch). 
This section contains the detailed performance results based on sequential 
processing of tracking data*. The error modeling approach and major find- 
ings were presented In Section 5. 

Cases for evaluation were defined originally for three TDAS constellation 
options (see Figure 5.2) and six user orbit types (see Figure 5.3), The 
results that follow cover 12 cases: four of the user orbits (A, D, E and 
F)** with each of the three TDAS constellations.** Each case Is presented 
In a standard format: 

• Tracking Schedule - This shows tracking contact intervals over 
24 hours for Schedules I and II and for beacon tracking (see 
Table 5-1). An accompanying user ground trace Indicates coverage 
outages due to a zone of exclusion (If any). 

• User Position & Time Accuracy Profiles - These show the la errors 
corresponding to the 24 hour tracking schedules and a subsequent 
6 or 12 hour prediction-only Interval. 

Summary tables (D-l thru D-4) Included at the end of Appendix D record 

the peak errors and Identify the major error contributors for each tracking 

alternative. 


•k ~ ~ 

Performance results for sliding batch processing are In Appendix E. 

** Individual cases are Identified by a two element mnemonic; e.g., A2 indi 
cates user orbit Type A and TDAS constellation Option 2. Evaluation of 
user orbit types, B (28°, 600 km) and C (28°, 1000 km) remains to be com 
pleted. 
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FIGURE D-3 

CASE A2 TRACKING SCHEDULES 
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APPENDIX E 


USER NAVIGATION PERFORMANCE RESULTS - SLIDING BATCH DATA PROCESSING 


TDAS user navigation performance was evaluated for two of the one-way track- 
ing alternatives (FLST and RLST) assuming sliding batch processing of track- 
ing data. The error modeling approach and major findings were discussed in 
Section 5; this section presents further details of the results obtained. 

Cases for evaluation were selected from among the TDAS constellation options 
and user orbit types defined in Figures 5-1 and 5-2. The results that fol- 
low cover half of the cases presented In Appendix D for sequential processing 
— the same user orbits types (A,D,E,F) and two of the TDAS satellite con- 
stellation options (1,2) . 

For each case considered, two types of user position and time accuracy pro- 
files were generated: 

• Single Batch Profiles - These were computed to assess the sensi- 
tivity of user QD/TD accuracy in the prediction mode to tracking 
interval selection. Durations of 6, 12, 18, and 24 hours were 
evaluated for each tracking schedule (1 and II).* Figures E-l 
and E-2 are examples. The tracking interval yielding the lowest 
peak error in the first 6 hours from end-of-tracking was selected 
for generating sliding batch profiles in each case. 


Schedules I and II for a 24 hour tracking interval are identical to those 
given in Appendix D for corresponding cases. Scheduled contacts for 
shorter tracking intervals were simply deleted. 


VI-E-1 


• Sliding Batch Profiles - These are comprised of prediction Inter- 
val segments corresponding to u sequence of single batch solutions, 
as defined earlier 1 ri Figure 5-4. Twenty-four hour profiles were 
constructed based on a batch Interval (P) of 3 hours and a com- 
putation/data handling Interval (C) of 1.5 hours. Figure E-3 
for Case A1 is an example.* 

Summary tables (E-l thru E-4) Included at the end of Appendix E record the 
peak errors and Identify the major error contributors for each tracking 
alternatl ve. 


The tlme-from-epoch scale on all plots has the same interpretation as In 
Appendix 0, l.e., the time elapsed since the user's initial location, 

{0°N , 0°E) . Individual batch epochs were adjusted backward or forward 
depending on tracking Interval In order to fit a given prediction segment 
slot. 
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FIGURE E-4 
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FIGURE E-ll 
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FIGURE E-19 

OSITI ON AND TIME ACCURACY WITH 1-WAY TDAS TRACKING - CASE F2 
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APPENDIX F 


GLOSSARY OF ACRONYMS 


BRTS 

Bilateration Ranging Transponder System 

CONUS 

Contiguous United States 

DSN 

Deep Space Network 

EOT 

End of Tracking 

FLBT 

Forward Link Beacon Tracking 

FLST 

Forward Link Scheduled Tracking 

GPS 

(NAVSTAR) Global Positioning System 

KSA 

K u -Band Single Access 

LSA 

Laser Single Access 

MA 

Multiple Access 

MCC 

Mission Control Center 

NASCOM 

NASA Cciiimuni cation Network 

NCC 

Network’ Control Center 

OD 

Orbit Determination 

OD/TD 

Orbit and Time Determination 

OSCF 

Orbit Support Computing Facility (GSFC) 

POCC 

Project Operations Control Center 

RLST 

Return Link Scheduled Tracking 

SA 

Single Access 

SMA 

S-Band Multiple Access 

TD 

Time Determination 

VLBI 

Very Long Baseline Interferometry 

WSA 

W-Band Single Access 

WSN 

White Sands (NM) 

A VLBI 

Differential VLBI (using two signal sources) 
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